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ABSTRACT 

We use the deep ground-based optical photometry of the Lyman Break Galaxy (LBG) Survey to 
derive robust measurements of the faint-end slope (a) of the UV luminosity function (LF) at redshifts 
1.9 < z < 3.4. Our sample includes > 2000 spectroscopic redshifts and w 31000 LBGs in 31 spatially- 
independent fields over a total area of 3261 arcmin 2 . These data allow us to select galaxies to 0.07L* 
and 0.10L* at z ~ 2 and z ~ 3, respectively. A maximum-likelihood analysis indicates steep values of 
a(z = 2) = —1.73 ± 0.07 and a(z = 3) = —1.73 ± 0.13. This result is robust to luminosity-dependent 
systematics in the Lya equivalent width and reddening distributions, is similar to the steep values 
advocated at z > 4, and implies that w 93% of the unobscured UV luminosity density at z ~ 2 — 3 
arises from sub-L* galaxies. With a realistic luminosity-dependent reddening distribution, faint to 
moderately luminous galaxies account for > 70% and > 25% of the bolometric luminosity density and 
present-day stellar mass density, respectively, when integrated over 1.9 < z < 3.4. We find a factor of 
8 — 9 increase in the star formation rate density between z ~ 6 and z ~ 2, due to both a brightening 
of L* and an increasing dust correction proceeding to lower redshifts. Combining the UV LF with 
stellar mass estimates suggests a relatively steep low mass slope of the stellar mass function at high 
rcdshift. The previously observed discrepancy between the integral of the star formation history and 
stellar mass density measurements at z ~ 2 may be reconciled by invoking a luminosity-dependent 
reddening correction to the star formation history combined with an accounting for the stellar mass 
contributed by UV-faint galaxies. The steep and relatively constant faint-end slope of the UV LF at 
z > 2 contrasts with the shallower slope inferred locally, suggesting that the evolution in the faint-end 
slope may be dictated simply by the availability of low mass halos capable of supporting star formation 
at z < 2. 
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1. INTRODUCTION 

The last decade has seen significant advances in our un- 
derstanding of the history of star formation and stellar 
mass assembly. Today, one can find several hundred de- 
terminations of the star formation rate density (SFRD) 
estimated from observations at many wavelengths across 
a large range of lookback time. Taken together, these ob- 
servations suggest a rapid increase in the SFRD from the 
epoch of reionization to z ~ 1 — 2, after which time the 
SFRD steadily decreased over the last ~ 10 billion years. 
This picture is generally understood in the context of hi- 
erarchical buildup at early times and gas exhaustion or 
heating at late times. While this characterization of the 
star formation history is broadly accepted, there are sev- 
eral key details that are missing from this interpretation, 
including the potentially important contribution of UV- 
faint (sub-L*) galaxies to the census of star formation 
and baryon budget. If rest-UV/optical light is a biased 
tracer of galaxy formation — particularly at high red- 
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shift where most of the baryons in galaxies are likely to 
reside in cold gas (jProchaska fc Tumlinsonll2"008h — then 
faint galaxies may constitute an important population for 
studying the process of star formation and feedback. Fur- 
ther, the number density of both bright and faint galax- 
ies departs significantly from expectations based on the 
ACDM model, suggesting a regulation of star formation 
in both luminosity regimes. In this paper, we present 
an extension of our previous work on the UV luminosity 
function at z ~ 2 — 3 in order to provide robust con- 
straints on the prevalence of UV-faint galaxies at epochs 
when galaxies were forming most of their stars. 

The luminosity function (LF) is a fundamental probe 
of galaxy formation and evolution, and can be used to 
address the relative importance of bright and faint galax- 
ies to the energy budget at a given epoch. Furthermore, 
comparison with the dark matter halo distribution in- 
forms us of the efficiency of star f ormation and effects of 
feedback at different mass scales (fRees k, OstrikerlH977t 
iSilk fc Reedll9M iDekel fc Birn boim 2001. Therefore, 
constraining accurately the shape of the LF is a necessary 
step in acquiring a more complete census of galaxies and 
elucidating the relationship between the baryonic pro- 
cesses that govern galaxy evolution and the dark matter 
halos that host them. 

The UV LF is relevant in several respects. Rest- 
frame UV emission is a direct tracer of massive star 
formation, modulo the effects of dust. Rest-frame UV 
observations of high redshift galaxies are generally not 
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limited by spatial resolution and the deepest observa- 
tions are up to a factor of ~ 2000 times more sensi- 
tive than those in the infrared and longer wavelengths. 
The combination of resolution, sensitivity, and the ac- 
cessibility of UV wavelengths over almost the entire age 
of the Universe makes the UV LF a unique tool in as- 
sessing the star formation history. The relative effi- 
ciency of UV-dropout selection has enabled a number 
of investigations of the UV LF at h igh redshifts based 
on photometrically-se l ected sa mples (ISteidel et alj | 1999 b 
Adelberger fc Steidell 120001 lYan fe Windhorstl 1200 



Bunker et al. 2004; Dic kinson et al. 2004; Bouwcns ~et all 
2007t ISawicki fc Thompson] 120061: [Reddv et alJ 120081 



Ouchi et alJ 12004 IBeckwith et all 120061: lYoshida et all 
20061: llwata et al.ll2007D . ~ 



Apart from the uncertainties that can be constrained 
from photometry alone, such as photometric errors and 
field-to-field variations, spectroscopy is a critical means 
of quantifying several important systematics that can af- 
fect the LF. These include contamination from low red- 
shift objects and high redshift AGN/QSOs, attenuation 
of UV emission due to the opacity of the intergalactic 
medium (IGM), and perturbation of galaxy colors due to 
Lya, reddening, and stellar population ages of galaxies. 
The relevance of these systematic effects is underscored 
by the fact that while there are numerous studies of the 
UV LF, the results have been inconsistent, both at low 
< 3; e.g.. lReddv et alj|2008t iLe Fevre ; et al.ll2005D and 
at hig h (z > 4: e.g.. IBeckwith et al . 2006: lBouwe~ et alJ 
120071 1 redshifts. 

Unfortunately, spectroscopic surveys have been limited 
to UV-bright (TZ < 25.5) galaxies at relatively low red- 
shifts (z < 4) and spectroscopic constraints on the num- 
ber density of faint sources at z > 2 are still lacking. 
Deeper spectroscopy is expensive and will remain so un- 
til the next generation of large (> 10 m) ground-based 
telescopes come online. Given the present practical lim- 
itations and the complexity of systematics involved in 
computing LFs, it seems prudent to revisit and extend 
our initial estimate of the UV LF (jReddv et all 120081 . 
hereafter R08) by evaluating the impact of these system- 
atics on the inferred number density of faint galaxies at 
z ~ 2 — 3. We go beyond the initial analysis of R08 
by quantifying several important effects relevant in the 
computation of star formation rate and stellar mass den- 
sities, including the effects of luminosity-dependent dust 
corrections and the integrated stellar mass of low mass 
galaxies. 

To this end, we combine what we know about the spec- 
troscopic properties of LBGs at z ~ 2 — 3 with deep 
ground-based optical data in 31 spatially independent 
fields to determine the faint-end slope with greater preci- 
sion. A brief description of the LBG survey, photometry, 
and followup spectroscopy is given in § [2] Our method 
for computing the UV LF is presented in § [3] and re- 
sults are discussed in § |4j Our results are compared with 
those in the literature and analyzed in the context of the 
early Hubble Deep Field (HDF) results, and we assess 
the impact of sample variance in § [5] The evolution of 
the UV LF is quantified in § [6] The contribution to the 
faint-end population from dusty, star-forming galaxies 
as well as quiescently-evolving galaxies with large stel- 
lar masses is discussed in § [JJ § [8] and [9] present con- 
straints on the star formation rate density (SFRD) and 



its evolution. We also reassess the stellar mass density 
at z ~ 2 and compare it to inferences from integrat- 
ing the star formation history. Finally, the evolution of 
the faint-end slope of the UV LF is discussed briefly in 
§ [10] All magnitudes are expressed in AB uni t s, un - 
less stated otherwise. Unless indicated, a lKroupal (|2001l ) 
IMF is assumed. A flat ACDM cosmology is assumed 
with H = 70 km s" 1 Mpc" 1 , fl A = 0.7, and Q, m = 0.3. 

2. DATA: SAMPLE SELECTION AND SPECTROSCOPY 

2.1. Fields and Photometry 

The Lyman Break Galaxy (LBG) survey is be- 
ing conducted in fields chosen primarily for hav- 
ing relatively bright background QSOs with which to 
study the interf ace between the IGM and galaxies 
at z ~ 2-3 (LAdelbgrger et all [20031 . 121)01120055 
ISteidel et afll2003L l2004h . Additionally, the survey was 
extended to include fields that are the focus of multi- 
wavelength campaigns, including the Groth- Westphal 
(|Steidel et alJ [2001: iGroth et al l [1991 and GOODS- 
North fields (Dickins on et al.1 120031 : iGiavalisco et al.1 
l2004bf) . Presently, the survey includes 31 fields, 29 of 
which have been targeted spectroscopically (see R08 for 
further details). As emphasized throughout this paper, 
the large number of spatially independent fields provides 
a precise handle on the magnitude of sample variance, 
an effect that has limited the conclusions that could 
be drawn from previous determinations of the LF. For 
this analysis, we have included two additional fields be- 
yond the 29 that were presented in R08, "Q1603" and 
"Q2240". In struments used and dates of observation are 
presented in ISteidel et al.l (|2003l[200l . For ease of refer- 
ence, the fields are listed in TableUJ together they include 
an area of 3261 arcmin 2 , or 0.91 square degrees. 

Excepting Q1603, a modified version of FOCAS 
(|Valdesl[l982l) was used to extract photometry from the 
optical ( U n GlZ) images of the su rvey fields. Source Ex- 
tractor (jBertin fc Arno uts 1996) was used for photom- 
etry in Q1603. We took care to minimize systematics 
between the FOCAS and Source Extractor results for 
Q1603 from an examination of the color distribution of 
recovered LBG candidates. Object detection was done at 
TZ band, and G — 1Z and U n — G colors were computed by 
applying the isophotal aperture from the TZ band image 
to the U n and G im ages. Further detail s on the photome- 
try ar e discussed in lSteidel et alJ (|2003h and lSteidel et al.1 
(|2004h . The images have a 5 a depth of 27.5 - 29.5 AB 
as measured through a ~ 3" diameter aperture. 

2.2. Color Selection 

We used the BX and LBG criteria which are based 
on the rest-frame UV colors expected of galaxies at red- 
shifts 1.9 < z < 2 . 7 and 2.7 < z < 3.4, respectively 
(|Steidel et al.ll2004 [20031 ). For spectroscopic followup, 
the sample of candidates was limited to TZ = 25.5. How- 
ever, because we are interested in using the entire pho- 
tometric sample to constrain the LF at 1.9 < z < 3.4, 
we did not impose this restriction. Rather, the detection 
significance and color distribution of candidates were ex- 
amined to establish a faint limit. The limits applied to 
each field and numbers of candidates are listed in Table 1. 
Fields with the deepest imaging in the survey allow us 
to select candidates to 1Z ~ 26.5 with a detection com- 
pleteness of w 60%, based on the simulations described 
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TABLE 1 
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TOTAL ... ... 3260.8 ... 22166 8663 



a Right ascension in hours, minutes, and seconds. 
b Declination in degrees, arcminutes, and arcseconds. 

c Number of BX candidates to the limiting magnitude given in column (5). 
d Number of LBG candidates to the limiting magnitude given in column (5). 



below. Our maximum-likelihood method for computing 
the LF allows us to extend the absolute magnitude limit 
w 0.5 mag fainter given the broadness of the redshift 
distribution, N(z), of the sample as constrained from 
extensive spectroscopic followup. Even at these depths, 
galaxies are detected with typically > 5 a significance 
in the 7^-band. The photometric sample used here in- 
cludes BX candidates in the original LBG survey fields 
where no spectroscopic followup of BX candidates was 
undertaken. 

2.3. Spectroscopy 

The spectroscopy of UV-selecte d candidates i nclud- 
ing LBGs and BXs is discussed in ISteidel et alJ (|2003l . 
l200l . To date, roughly 24% and 35% of BX and LBG 
candidates, respectively, with IZ < 25.5 have been tar- 
geted spectroscopically. The resulting sample includes 
2023 star- forming galaxies with 1.9 < z S pcc < 3.4, the 
largest spectroscopic sample of star-forming galaxies at 
these redshifts. The spectroscopic statistics, including 
the number of spectroscopic redshifts, for each field of 
the LBG survey are given in R08. 

The spectroscopic sample is used to estimate the over- 
all contamination rate in the photometric sample. This 
contamination can arise from stars, low redshift galaxies 
and AGN, as well as QSOs and AGN at 1.9 < z < 3.4. 
The contamination statistics and the extent to which 



they can be applied to determine the number of contam- 
inants in the photometric sample are discussed in R08. 
Note that the contamination rate is a strong function 
of magnitude (being largest at bright magnitudes) and 
quantifying it is a crucial step in computing the bright- 
end of the LF. 

3. INCOMPLETENESS CORRECTIONS 
3.1. Method 

The approach that we adopted to correct the LBG 
sample for incompleteness and compute the UV LF is 
described in detail by R08. For convenience, here we 
summarize a few of the key features of our method. The 
primary goal is to construct a set of transformations that 
relate the observed properties of galaxies (e.g., observed 
luminosity, rest-frame UV slope, and redshift) to their 
true properties (e.g., intrinsic luminosity, reddening, and 
redshift). Using X-ra y and mid-IR data fo r a sample of 
LBGs at z ~ 2 - 3, iReddv et all (|2006bl ) demonstrate 
that the rest-frame UV slope can be used to measure the 
amount of dust reddening for typical LBGs, and we will 
assume this for the subsequent discussion. 

We first used a Monte Carlo simulation to add galax- 
ies of varying sizes and colors to our U n GTZ images. 
The distribution of colors reflects that expected for star- 
forming galaxies at z ~ 2—3 with constant star formation 
for > 100 Myr and varying amounts of dust reddening. 
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Specifically, we added objects with redshifts 1 < z < 4 
and reddening of 0.0 < E(B - V) < 0.6. The simulated 
rcdshift i s used to appl y an IGM opacity to the colors 
using the Madau (1995) prescription. To make the sim- 
ulation as realistic as possible, we force d the simulated 
galaxies to abide by a ISchechterl (|1976l ) luminosity dis- 
tribution and added just 100 — 200 of them at a time 
to the images. The latter restriction maintains the de- 
blending statistics which in turn affect the photometric 
errors. Photometry was performed in the same man- 
ner as was used on the real data and the detection rate 
and recovered magnitudes and colors of the simulated 
galaxies were recorded. We repeated this procedure un- 
til k 2 x 10 5 galaxies were added to the images in each 
field of the survey. 

It is common to use such simulations to determine what 
fraction of galaxies with a given magnitude will be de- 
tected with colors that satisfy the LBG selection crite- 
ria. However, strictly speaking, the simulation will only 
tell us the probability that galaxies with a given simu- 
lated magnitude will be detected with a given recovered 
magnitude, and there may not be a monotonic corre- 
spondence between simulated and recovered magnitude. 
More generally, it is not necessarily true that the av- 
erage simulated properties of galaxies arc equivalent to 
their average observed properties. This is particularly 
true if photometric errors have significant biases and are 
comparable to the bin sizes used to compute the LF and 
the selection window spans a region of color space not 
much larger than the typical photometric errors. Other 
systematic effects, such as the Lya equivalent width dis- 
tribution (Whya) of the population, may scatter galaxies 
in certain directions of color space. Also, some galax- 
ies with simulated colors that do not initially satisfy the 
color criteria may have recovered colors that do: by def- 
inition, these galaxies' simulated properties will not, on 
average, reflect their observed properties. Because of 
these systematic effects, the number of galaxies that lie 
in a particular bin of observed properties will be some 
weighted combination of the number of galaxies in any 
number of bins of simulated properties. 

Because of these systematic effects, we approach the 
problem of incompleteness by using the maximum like- 
lihood method described in R08. Using this formalism, 
our goal is to maximize the likelihood of a given set of lu- 
minosity, reddening, and redshift distributions, denoted 
by C, according to the following expression: 
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where n^k is the expected number of galaxies in the 
i th bin of luminosity, j th bin of reddening, and k th bin 
of redshift that the values of the luminosity, reddening, 
and redshift distributions imply and mjk is the observed 
number of galaxies in that bin. The Monte Carlo sim- 
ulations give the probability that a galaxy in the i'j'k' 
bin of simulated properties will lie in the ijk bin of re- 
covered properties. The set of probabilities, defined as 
the transitional probability function, 
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is used to compute n^k- The basic procedure is to then 
vary the input distributions of luminosity, reddening, and 



Fig. 1. — Intrinsic rest-frame Lya equivalent width (ll^Lya) dis- 
tribution for 1Z < 25.5 star-forming galaxies at 1.9 < z < 2.7 
and 2.7 < z < 3.4, from R08. Lya in emission is represented as 
Wi^ya > 0. Error bars are determined from simulations and reflect 
the variance in the distributions allowed by the errors in the UV 
LF and reddening distribution (see R08 for discussion). 

redshift until the differences between the expected and 
observed numbers of galaxies in each ijk bin are mini- 
mized. 

3.2. Lya Equivalent Width (Wi^ ya ) Distribution 

An important systematic effect to consider is the scat- 
tering of colors due to the presence of Lya emission and 
absorption: the Lya line falls within the U n and G bands 
at redshifts that are targeted with the BX and LBG cri- 
teria, the same bands that are used to select the galaxies. 
Rather than adding the W^ya distribution as another free 
parameter in the maximum-likelihood analysis — thus 
complicating our ability to determine the LF — we per- 
formed simulations where we made various assumptions 
regarding the intrinsic W^ ya distribution at z ~ 2 — 3 and 
observed the effects on the best-fit LF (R08). We did this 
by adding a random W\^ ya drawn from a distribution of 
W~L ya to each galaxy, and then recomputing the colors of 
each galaxy. Effectively, the addition of Lya will perturb 
the colors and thus modulate the transitional probability 
function, £. R08 showed that the BX and LBG color cri- 
teria did little to alter the intrinsic Wh ya distribution at 
z ~ 2 and z ~ 3. Therefore, we assume the W\^ ya distri- 
bution observed for BXs and LBGs (Figure Q}. Here, we 
repeat the simulations of R08, but also allow for changes 
in the shape of the Wz, ya distribution proceeding from 
UV-bright to UV- faint galaxies (see Appendix). 

4. RESULTS: THE UV LF AT 1.9 < Z < 3.4 

4.1. Computation of the LF and Errors 

The value of the luminosity distribution that maxi- 
mizes the likelihood of observing our data (Eq. [1]) is com- 
puted using the method discussed in the previous sec- 
tion. Initially, we assumed that the intrinsic distribution 
of (1) rest-frame UV slopes, (2) redshifts, and (3) W\, ya 
remain constant as a function of apparent magnitude. 
In R08, we justified these assumptions when computing 
the LF to our spectroscopic limit of 1Z — 25.5. How- 
ever, it is not unreasonable to suspect that, for example, 
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the reddening and W\^ ya distributions of galaxies fainter 
than our spectroscopic limit may be different than those 
for galaxies where we are able to directly measure the 
distributions with spectroscopy ( "spectroscopic distribu- 
tions" ). Such differences will change £ and thus affect our 
incompleteness corrections. First, we first made the sim- 
plified assumptions that all of these distributions remain 
unchanged as a function of apparent magnitude. LFs de- 
rived in this case are referred to as the "fiducial" LFs. 
In the appendix, we discuss how the LF would change 
if we adopt more realistic assumptions for the properties 
of UV-faint galaxies. In our analysis, the effect of in- 
creasing photometric error for fainter galaxies is already 
incorporated in the calculation of £. 

LFs were computed separately for star-forming galax- 
ies in the redshift ranges 1.9 < z < 2.7 and 2.7 < z < 3.4 
using the photometric BX and LBG samples, respec- 
tively. For the lower redshift range, the LF is computed 
in terms of a (composite) absolute magnitude that is the 
average of the G and 7?.-band fluxes. For the higher red- 
shift range, the LF is computed using the 7?.-band mag- 
nitude. This method provides the closest match between 
rest-frame wavelengths, roug hly 1700 A. 

The total error in the LF is computed using the follow- 
ing method. The observed number counts of galaxies in 
each field were adjusted randomly in accordance with a 
Poisson distribution and the maximum-likelihood LF was 
computed for each field. This procedure was repeated 
many times for all the fields. The dispersion in the LF 
values for each bin in absolute magnitude is taken as the 
total error which, as a consequence of the way in which 
it is computed, includes both Poisson and field-to- field 
variations. 

4.2. Summary of Systematic Effects and Final Results 

The details of the systematic tests performed to judge 
the effects of luminosity-dependent W^ ya and reddening 
distributions on the LF are presented in the Appendix. 
To summarize, we analyzed the influence of galaxies with 
(1) strong Lya emission and (2) zero or declining redden- 
ing with apparent UV magnitude. Employing current 
estimates of the mean stellar population, reddening, and 
number density of galaxies with strong Lya emission as a 
function of UV luminosity at high redshifts, we find that 
the inferred number density of galaxies on the faint-end 
of the UV LF increases by < 3% at 1.9 < z < 2.7 and de- 
creases by < 4% at 2.7 < z < 3.4. Because these changes 
are not negligible compared to the Poisson and field- 
to-field errors on the faint-end number densities, they 
should be included in any proper assessment of the LF. 
Nonetheless, these changes in number density can be ac- 
commodated by Schechter parameterizations that vary 
within the uncertainties of the individual parameters, a, 
M*, and <p*. 

We have also examined how changes in the mean red- 
dening of galaxies as function of UV luminosity affects 
our measurement of the UV LF. We considered two sce- 
narios, one in which the extinction drops to zero for 
galaxies fainter than 1Z — 25.5 and one in which the 
extinction decreases monotonically with UV luminosity, 
and approaches zero in the faintest luminosity bin con- 
sidered here. The latter scenario is more realistic than 
the former, and is parameterized as a linear relation be- 
tween E(B — V) and magnitude (see the appendix; we 



TABLE 2 

Rest-Frame UV Luminosity Functions of 1.9 <z < 3.4 Galaxies 



4> 

Redshift Range M A b(1700A) (xlO -3 /ig 7 Mpc -3 mag -1 ) 



-22.83 — 


-22.33 


0.004 ± 0.003 


-22.33 — 


-21.83 


0.035 ±0.007 


-21.83 - 


-21.33 


0.142 ±0.016 


-21.33 - 


-20.83 


0.341 ±0.058 


-20.83 - 


-20.33 


1.246 ±0.083 


-20.33 - 


-19.83 


2.030 ±0.196 


-19.83 - 


-19.33 


3.583 ±0.319 


-19.33 - 


-18.83 


7.171 ±0.552 


-18.83 - 


-18.33 


8.188 ±0.777 


-18.33 - 


-17.83 


12.62 ± 1.778 


-23.02 — 


-22.52 


0.003 ± 0.001 


-22.52 — 


-22.02 


0.030 ±0.013 


-22.02 — 


-21.52 


0.085 ± 0.032 


-21.52 — 


-21.02 


0.240 ±0.104 


-21.02 — 


-20.52 


0.686 ± 0.249 


-20.52 — 


-20.02 


1.530 ±0.273 


-20.02 - 


-19.52 


2.934 ± 0.333 


-19.52 - 


-19.02 


4.296 ± 0.432 


-19.02 - 


-18.52 


5.536 ±0.601 



TABLE 3 

Best-fit Schechter Parameters for UV LFs of 1.9 < z < 3.4 
Galaxies 



Redshift Range a M^ B (1700A) <j>* (xlO -3 Mpc -3 ) 

1.9 < 2 < 2.7 -1.73 ±0.07 -20.70 ±0.11 2.75 ±0.54 

2.7<z<3.4 -1.73 ±0.13 -20.97 ±0.14 1.71 ±0.53 



refer to this latter scenario as the "luminosity-dependent 
reddening model"). In this case, we find appreciable 
increases of « 10% in the inferred number density be- 
tween 1.9 < z < 2.7. In the higher redshift range 
2.7 < z < 3.4, there is little change in the inferred num- 
ber densities. Our final determinations of the LF and 
the corresponding Schechter fits are shown by the data 
points and solid lines, respectively, in Figure [U with the 
values and Schechter parameterization listed in Tables [2] 
and[3l 

Our determinations of the bright-end of the UV LFs 
to M (17001) = -18.83 and -19.52 at z ~ 2 and z ~ 3, 
respectively, incorporate data over 3261 arcmin 2 in 31 
independent fields. Data from 22 fields and 2239 arcmin 2 
are used to constrain the LFs at -18.83 < M(1700A) < 
-18.33 and -19.52 < M(1700l) < -19.02 at z ~ 2 
and z ~ 3, respectively. Finally, data from 6 spatially 
independent fields and 317 arcmin 2 are used to constrain 
the LF in the faintest magnitude bin to M(1700A) = 
— 17.83 and —18.52 at z ~ 2 and z ~ 3, respectively. To 
ensure that spatial variance in these 6 deep fields are not 
driving the observed faint-end slope, we recalculated a 
by fitting Schechter functions to the LFs excluding the 
faintest bin. Allowing <f>* and M* to vary, we calculate 
a = -1.75 ± 0.09 and -1.94 ± 0.18 at z ~ 2 and z ~ 3, 
respectively, still significantly steeper than the shallower 
a > —1.6 found in previous studies. The similarity in a 
obtained with or without data from the faintest bin is not 
surprising given that the uncertainty in the LF includes 
the sample variance from the 6 fields used to constrain 
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Fig. 2. — Rest- frame UV luminosity functions of star-forming galaxies at 1.9 < z < 2.7 (circles) and 2.7 < z < 3.4 (squares), along with 
the best-fit Schcchter (1976) functions. The 68% and 95% likelihood contours between M* and a for our final determinations of the LFs 
are shown in the inset panel. 

the number density in this bin. 

The degeneracy between the faint-end slope (a) and 
characteristic magnitude (M*) — illustrated by the like- 
lihood contours in Figure [5] (inset) — is reduced signifi- 
cantly compared to that computed in R08. Our analysis 
extends to luminosities that are 4 times fainter than the 
limit dictated by efficient spectroscopy, and « 14 and 10 
times fainter, respectively, than the characteristic lumi- 
nosity L* at z ~ 2 and 3. Our sample is large enough so 
that the error in the LF at all magnitudes is dominated 
by field-to-field variations (Figure [3|) . Within the total 
errors, the UV LFs at z ~ 2 and z ~ 3 are virtually in- 
distinguishable, indicating little change between the two 
in the number density of both UV-bright and UV-faint 
galaxies. 

5. DISCUSSION: LARGE-SCALE CONTEXT 

In the following sections, we discuss our results in 
the context of previous determinations of the UV LF 
(§ 15. ip and its evolution with redshift (§ [6j) . To gain 
further insight into the nature of sub-L* galaxies, we as- 
sess the contribution to the faint-end population from 
dusty, star-forming galaxies and those with large stel- 
lar masses (§0. In § [51 we discuss the implications of 
a luminosity-dependent reddening distribution and the 
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Fig. 3. — Poisson (dashed line) and total (solid line) error in the 
UV LF at z ~ 2. The Poisson error increases to fainter magnitudes 
given the smaller survey area used to constrain the faint-end. At 
all magnitudes, however, the error in the LF is dominated by ficld- 
to-ficld variations. Similar results are obtained for the z ~ 3 UV 
LF. 

average corrections required to recover the bolometric 
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Fig. 4. — Comparison o f UV LFs at z ~ 2 (left) a nd z ~ 3 (right). For clarity, data points are excluded on all but our current 
determinations and those of Sawicki & Thompson (2006), but the errors are typically smaller than the observed differences discussed in the 
text. Also shown are the survey areas over which the LF is derived, with some surveys using a combination of wider shall ower data to ancho r 
the bright-end of the L F and deeper data in smalle r areas to constrain th e faint-end slope. Includ e d are are data fromlReddv et all (20081: 
ISawicki fc Thompson! i200^) ; IPaltani et all ll2007h ; lArnouts et"aT] ll2"005h ; IGabasch et~aTI (12001) ; IPoli et al.l l EfflOlI) ; lAdelberger fc Steidell 
[pOOOMSteidel et alj 1)19991) . 

star formation rate density. These findings are then dis- 
cussed in the context of the star formation history and 
the buildup of stellar mass in § Finally, we discuss 
briefly the evolution of the faint-end slope in § [TO] 



5.1. Comparisons at z ~ 2 — 3 

In this section, we place our results in the context of 
previous determinations of the UV LFs, starting with 
those around z ~ 2. Figure [4] summarize s the results of 
severa l previous studies including th ose of Gabasc h et al.l 
(|2004l ); ISawicki fc Thompson! (|2006l ) and R08, along with 
our current determination, at z ~ 2. The redshift inter- 
vals over which the LF is computed are similar between 
these studies, but we note the almost two orders of mag- 
nitude difference in the areas probed, from rj 40 arcmin 2 
at the low end to 3261 arcmin 2 in the current determina- 
tion. There are significant differences between the LFs 
at faint magnitudes (M(1700A) > -20). In general, it 
is possible that the determinations of t he smaller sur- 
veys (e.g., from the FORS-Deep Field; IGabasch et al.l 
2004) could be mimicked by an overall under-density 
in the small area probed combined with an overestima- 
tion of the bright-e nd due to contamina tion from low 
redshift interlopers. IGabasch et al.l (|2004h do not spec- 
ify the contamination fractions for their higher redshift 
samples at z > 2, so a fair comparison with our find- 
ings at the bright-end is not possible. We also note 
that IGabasch et alj (|2004l ) relied on photometric red- 
shifts which could not be well-calibrated due to the lack 
of spectroscopically-confirmed galaxies at z ~ 2 exam- 
ined in their study (see their Figure 2). R08 showed 
that biases in photometric redshifts can easily boost the 
bright-end of the luminosity function with respect to the 
faint-end, and the overall shape of their LF may be a 
result of this effect. 

Perhaps a fairer compa rison can be made with 
ISawicki fc Thompson! (2006) since they use exactly the 
same filter set to select BX candidates at z ~ 2.3 in the 
Keck Deep Fields (KDFs). For their fiducial model, they 
assumed no perturbation of colors due to Lya and a con- 
stant E(B — V) = 0.15 with no dispersion, and they com- 
pute the LF using the standard V e s method. Their LF 
suggests a much shallower slope of a ~ —1.2 compared 



to our result (Figured]). However. ISawicki fc Thompson! 
(2006) point out that their LF derived at z ~ 2 is sen- 
sitive to the assumed E(B — V), and that bluer val- 
ues of E(B — V) will tend to yield larger inferred num- 
ber densities (see their Figure 7). This observation is 
consistent with our finding and, in particular, if the 
E(B — V) distribution becomes significantly bluer pro- 
ceeding to fainter galaxies, this effect would manifest it- 
self as a steepening of the faint-end slope (appendix). 
However, this systematic effect alone cannot account 
for all of the differen c e betw een our result and that of 
ISawicki fc Thompson! (2006), since even in the fiducial 
case of a luminosity-invariant E(B — V) distribution (but 
not constant- valued) we find a steep a — —1.67 ± 0.06 
(appendix). In any case, regardless of how the mean 
E(B — V) varies with magnitude, the distribution itself 
is not a delta function, of course, and has intrinsic dis- 
persion; those galaxies at the blue end of the distribution 
(i.e., less reddening) will tend to escape the selection cri- 
teria more frequently than galaxies with redder E{B—V) . 
Hence, a bluer mean E{B — V), the intrinsic scatter in 
E(B — V) for UV-faint galaxies, and a general perturba- 
tion of colors due to Lya will all result in larger corrected 
number densities at the faint-end. 

Focusing on the higher redshift range, we find again 
reasonable agreement among the various determinations 
of the bright-end of the UV LF at z ~ 3. 5 The 
only significantly discrepant points are from the VVDS 
that imply signific ant numbers of UV-bright galaxies 
(jPaltani et alJ l2007). However, applying the correct con- 
tamination fractions (based on spectroscopy) to their 
points brings them in accordance with the other deter- 
minations (R08). As at lower redshifts, we find a sub- 
stantially steeper faint-end slope at z ~ 3 than sug- 
gested by previous results. Most determinations have 
foun d a > —1.6, shallow er than the canonical a = —1.6 
from lSteidel et all (Il999h, alth ough most of these studies 
(including ISteidel et al.l Il999h constrained a using deep 
data from only one or two small deep fields (e.g., Hubble 
Deep Field, FORS Deep Field) where large-scale struc- 



° Sawicki & Thompson (2006) use the results from Steidel et al. 
(119991 ) to constrain the bright-end of the UV LF at z ~ 3. 
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ture may be an issue. ISawicki fc Thompson! (|2006l ) find 
a = —1.43^09 based on the Keck Deep Field data over 
an area of 169 arcmin 2 . 

5.2. Differences in LF Computation 

What could be the reason for the disparity in the 
faint-end number densities between our study and pre- 
vious determinations? Without a more detailed com- 
parative analysis incorporating the data used in these 
other studies, it is difficult to pinpoint a single cause 
for the discrepancy. There are, however, a number of 
differences between our analysis and others that may 
lead to the observed variance in a. Our analysis (1) 
uses over 2000 spectroscopic redshifts to evaluate and 
correct for contamination as a function of luminosity; 
(2) models the systematic effects of a luminosity depen- 
dence in the intrinsic Lya equivalent width and redden- 
ing distribution of galaxies, likely to be the two domi- 
nant sources of systematic error in the LF; (3) employs 
a maximum-likelihood method that is more robust than 
the V e g method against biases in photometry and other 
non-uniform sources of scatter; and (4) takes advantage 
of data in 31 spatially uncorrelated fields over a total area 
of close to a square degree. Even at the faint-end, our 
determinations are based on 6 independent fields with a 
total area of 317 arcmin 2 , a roughly 88% larger area than 
used in the previous faint-end determination at z ~ 2 — 3 
(but see next section) . For all of these reasons, we believe 
our LFs to be the most robust determinations to date. 

The differences in faint-end slope derived between 
studies with similar depth is not particularly significant 
within the margina lized errors on a. For exam ple, the 
a = -1.43+g;^ of ISawicki fc Thompson! (l2006h is still 
consistent within the 1 er (marginalized) error of our de- 
termination of a — —1.73 ± 0.13 at z ~ 3. Yet, the 
difference in the actual number density of faint galax- 
ies is significant at the 2 — 3 a level. This emphasizes 
why comparisons between a derived in different stud- 
ies should perhaps not be taken too seriously without 
placing them in the context of the errors on the actual 
number density of UV-faint galaxies. 

5.3. Cosmic Variance 

In spite of the care used in the present sample, even 
317 arcmin 2 is a relatively small area over which to con- 
strain a. As noted above, the uncertainties in the LF are 
dominated by field-to-field variance at all magnitudes. 
We can assess how the empirically-constrained errors on 
the UV LF compare to expectations based on the cor- 
relation function. Following the procedure outlined by 
iTrenti fc Stiave lli (2008!), we can estimate the combined 
uncertainty due to cosmic variance and Poisson statis- 
tics by integrating the two point correlation function for 
dark matter halos with some average galaxy bias. The 
basic premise is that the spatial correlation function of 
halos gives information on the variance in the spatial dis- 
tribution of galaxies along different lines of sight given 
various assumptions for the cosmology and halo filling 
factor. For this calculation, we assumed a number den- 
sity of objects as implied by the maximum-likelihood LF 
at z ~ 2 and a sample "completeness" fraction of 0.47. 
This number is the ratio of star-forming galaxies that 
satisfy the color selection criteria to the total number of 



star-forming galaxies as determined from the LF (see R08 
for a discussion of this fraction) . The cosmology is set as 
follows: VL X = 0.74, O m = 0.26, H Q = 70 k m s" 1 Mpc" 1 
a spe ctral index n s = 1, and a% — 0.9 (jSpergel et al.l 
|2007| ). We must also make some assumption for the 
halo filling factor. Star-forming galaxies are scattered 
out of the LBG selection window due to p rimarily ran- 
dom processes such as photometric scatter (|Reddv et al.l 
2005) , and must therefore cluster in the same way as 
galax ies that do satisfy the LBG criteria (|Conrov et al.l 
2008). Further, the comoving number density of LBGs 
is similar to the number density of halos that have simi- 
lar cl ustering strength, s uggesting a halo filling factor of 
re 1 (|Conrov et al.l l2008h . Assuming this remains valid 
for UV-faint galaxies, we find a fractional error in num- 
ber counts of re 9% over a survey area of 317 arcmin 2 . 
In general, we would expect this calculation to yield a 
lower limit to the uncertainty since other effects (e.g., 
uncertainties in zeropoints and systematics in the color 
distributions from field to field) contribute to the error 
in the LF, and indeed our empirically-derived error in 
the faintest bin of the z ~ 2 UV LF is re 50% larger 
than the value obtained from the two point correlation 
function. For comparison, this calculation implies that 
the field-to-field variance is re 17% lower than what we 
would have obta ined over the area probed b y the KDFs 
of 169 arcmin (|Sawicki fc Thompsodl200l . This dif- 
ference is not large enough to explain the apparent dis- 
crepancy at the faint-end, and some of the systematics 
discussed above are also likely to play a role. Deep UV 
imaging over areas of close to a square degree (similar 
to that used to estimate the bright-end of the LF) will 
be required to constrain the fractional error in number 
counts to < 5%. 

5.4. Hubble Deep Field (HDF) 

A comparison between the present work and those 
of the early HDF-based studies of the UV LF is use- 
ful, particularly in light of the often-used argument 
that the HDF presents a biased view of the Uni- 
verse, and one that is invoked to expla in the diver- 
gent r esults on the UV LF at z ~ 3 (ISteidel et alJ 
19991 iDickinson et all 120031: iGiavalisco et al.l 12004a; 
Gabasch et al.1 12004 ISawicki fc Thompson! 1200(1. The 
bright -end of the LFs computed here and bv Stei del et al.l 
((1991 are in excellent agreement. Within the 1 a 
marginalized errors, the faint-e nd slope derived at z ~ 3 
agrees with the slope found by ISteidel et al.l (|1999l ) and 
lAdelberger fc Steidell (|2000j ) . and there is essentially no 
significant difference i n 4>* and M*. Howe ver, given the 
widespread use of the ISteidel et al.l (|1999j ) results, it is 
important to note that their determination of a — con- 
strained from a [/-dropout sample in the HDF — does 
not take into account incompleteness from photometric 
scatter. As discussed in § G2 the effect of such scatter 
is to make the incompleteness corrections larger at the 
faint-end, thus steepening the faint-end slope. In sum- 
mary, contrary to the suggestion that the HDF contained 
an over-density of faint galaxies relative to bright ones 
when compared with other fields, our results imply that 
the HDF is reasonably representative of the z ~ 2 — 3 
universe. 

6. DISCUSSION: EVOLUTION OF THE UV LF 
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Fig. 5. — (Left) Evolution of the UV LFs from z ~ 7 to z ~ 2. For clarity and consistency, we show only LFs at z > 4 from Bouwcns et al. 
(t2007l . 2008) since they are calculated using a maximum-likelihood technique similar to the one used here. For comparison, the local UV 
LF from Wydcr ct al. (2005) is also shown. (Right) E volution of the charac teristic UV luminosity or magnitu de, M* , with r edsh ift. Points 
are from lWvder et al.l (120051 ) at z ~ (open triangle), Arnouts et al. (2005) at < z < 3.0 (filled tr ianales) , IBouwens et al.l (120081) at z > 4 
(squares), and our determinations at z ~ 2 — 3 (circles). 
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Fig . 6. — (Left) dN/dz as a function of redshift, assuming our determinations of the UV LF at z ~ 2 — 3 and those of (Bouwcns ct al. 
|2007f ) at z > 4, integrated to 0.1L*_ 3 (blue) and 0.1L* (red). The shaded regions indicate approximately the uncertainty based on the 
errors in the Schechter parameters. (Right) Total dL/dz as a function of redshift (green) and dL/dz brighter and fainter than 0.1L*_ 3 (red 
and blue, respectively). 



Figure [5] summarizes our UV LFs at z ~ 2 — 3 
along with higher redshift determinations. For clar- 
ity and consis t ency, we included the findings from 
IBouwens et al.l (|2007f) only since those authors use a 
maximum-likelihood method for determining the LF that 
is similar to the method we use. These authors provide 
a detailed comparison of UV LFs at z > 4 from different 
studies. 

6.1. Evolution in M* 

Despite the large number of investigations at z > 4, 
there is still a fair amount of uncertainty regarding the 
parameterization of the evolution in the UV LF. Some 
have claimed that the evolution occurs primarily in L* 
(IBouwens et al.l I2007D. while others find an evolution 
in </>* (|Beckwith et al.l 120061 ) or the faint-end slope a 
(|Iwata et al.l 12007^ Some have also suggested an evo- 
lution in both L* and </>* (e.g., iDickinson et all 120041 
iGiavalisco et al.l 1200 4a) such that the total luminosity 
density could remain constant from z ~ 3 — 6. Of course, 
the reach of some of these conclusions is limited by the 
depth of data used to derive the LF. Because the LFs 
at z > 2 shown in Figure [5] are derived using data of 
comparable depth and analyzed in a similar manner — 



although we note that our LFs at z ~ 2 — 3 are anchored 
by spectroscopy and photometry in an area roughly an 
order of magnitude larger than used at z > 4 — here- 
after we will assume that the evolution of the LF at 
z > 4 can be accommodat e d by a change in L* as ad- 
vocated by IBouwens et al.l (|2007l ) . Given the observed 
fading of galaxies at z < 2 (e . g.. [Dickinson et al.l 120031; 



lading ot gal axies at z >, 2 l e - g -- ILh ckmson ct al. UUUdt 
Madau et al]ll996l : iLillv et al.l fl995. 1996; StcidcLet_alJ 



1999), it is useful to examine our results in the context of 



this evolution in L* (Figure[5]). In particular, we find that 
L* is brightest at z ~ 2 — 3, with this average unobscured 
UV luminosity decreasing at z > 4 (earlier cosmic time) 
and decreasing by a factor of ~ 16 between z ~ 2 and 
the present-day. Quantitatively, IBouwens et al.l (|2008l ) 
found a linear parameterization between M* and z at 
z > 4 that appears to follow that generally expected for 
the growth of the halo mass function — assuming an evo- 
lution in the mass-to-light ratio for halos of ~ (1+z) -1 — 
given standard assumptions for the matter power spec- 
trum, indicating that hierarchical assembly of halos may 
be dominating the evolution in M* , or equivalently L* . 
In the context of this study, the linear parameterization 
can be ruled out at the 8 a level at z — 2.3 (in the 
sense that it would predict a significantly larger L* at 
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Fig. 7. — Faint-end slope a as a functi on of redshift. At z > 2, we 
include our results ( filled circles), th at of Steidel ct al. ( 1999) (open 
circle), and those of Bouwcns ct al. (2007) (squares). At lower red- 
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shifts, we included only points in the IRvan et al.l 1120071) compila- 
tion that were derived from the rest-UV LF and that relied on 
data extendi ng at least two magnitu des f ainter than M* , incl uding 
results from Budavari et al. (2005|) and IWvder et all (120051) (tr i- 
angles). Also show n are points from lSawicki &: T hompson (2006), 
llwata et al.l (120031) (errors in a are no t provided by these authors; 
crosses), and Yan & Windhorst (2004) (range of likely a indicated 
by hashed box). The dashed line marks the mean value of a 
found at z > 2 from our study and that of Bouwcns ct al. (2007) 
((a) ~ -1.73). 

z = 2.3 than is observed), indicating that by these red- 
shifts, some other effect(s) modulate L* away from the 
value expected from pure hierarchical assembly. 

These observations are illustrated more clearly by ex- 
amining dN/dz as a function of redshift (Figure^, which 
is extrapolated based upon linearly fitting the relation- 
ship between L* and z and <jf and z, and assuming a 
fixed a = —1.73 as indicated by the Schechter fits (Ta- 
ble [3]). Integrating the number counts to a fixed luminos- 
ity shows that bright galaxies with L > 0.1L* =3 increase 
in number density by an order of magnitude with cosmic 
time from z ~ 7 to z ~ 2. Alternatively, the number 
counts are flatter at z > 4 when integrating to Q.\L*{z) 
(i.e., L* appropriate at the redshift z where dN/dz is 
calculated) suggestin g that (j>* is relativel y constant at 
these redshifts (e.g., iBouwens et al.l [20071) . There may 
be a slight increase in 4>* between z ~ 2 — 3, though the 
magnitude of the errors on </>* are large enough that we 
cannot rule out non-evolution in the number density. 

Also shown is dL/dz, both above and below a fixed 
luminosity, in this case 0.1L* =3 , along with the total lu- 
minosity density. The evolution implied by our LFs sug- 
gests that the approximately order of magnitude increase 
in luminosity density between z ~ 7 and z ~ 4 is followed 
by a flattening betwe en z ~ 2 — 3. This res ult itself is 
hardly surprising (see lGiavalisco et "aT1l2004al R08), but 
its significance is constrained robustly given that our LFs 
are determined over two orders of magnitude in luminos- 
ity. We will return to a discussion of these findings in 
the context of the cosmic star formation history (§|SJ). 

6.2. Evolution in a 

Perhaps the most striking result of our analysis — and 
one that is possible to address with confidence given the 
depth of data considered here — is a very steep faint- 
end slope of a ~ —1.73 at z ~ 2 — 3 that is robust to the 
luminosity-dependent systematics discussed in the Ap- 
pendix. The a we derive at z — 2.30 is virtually identical 
to that derived at z = 3.05, and is remarkably similar to 



the steep faint-end slopes favored at z > 4 (Figure [7]). 
Given the rapid evolution in L* and the luminosity den- 
sity at z > 2, the invariance of a over the s ame ~ 3 Gyr 
timespan and the shallow a found locally (|Wvder et al.1 
120051 : Budavar i et al.ll2005f ) pose interesting constraints 
on models of galaxy formation. We revisit this issue in 

§HJ3 



7. 



DISCUSSION: NATURE OF GALAXIES ON THE 
FAINT-END OF THE UV LF 



Before proceeding to discuss the implications of our 
results, it is useful to assess the contribution of galax- 
ies selected with different methods to the UV LF. R08 
demonstrate that the BX and LBG criteria to 1Z = 25.5 
select the majority of galaxies on the bright-end of the 
UV LF, namely those with Ltjv ^ 0.1L*. We show in 
the appendix that these criteria recover the majority of 
star- forming galaxies fainter than 0.1L*. The tests dis- 
cussed in the Appendix assume that the vast majority 
of galaxies on the faint-end of the UV LF are relatively 
unreddened, young galaxies. The aim of this section is to 
quantify the fraction of galaxies on the faint-end that are 

(1) UV- faint simply because they are heavily- reddened or 

(2) older galaxies that have passed their major phase of 
star formation. The latter investigation is relevant if we 
are to make inferences on the connection between the 
dark matter halo mass distribution and the luminosity 
function. 

7.1. Bolometrically- Luminous Galaxies 

Deep mid-to-far IR surveys have uncovered a siz- 
able population of dusty and i nfrared luminous 

galaxies at z ~ 2 — 3 (e.g. , lYan et ail 120071; 

Caputi et ail 120071: [Papovich et alj 120071: iReddv et all 



20051: IChapman et all 120051: Ivan Dokkum et al.l 12004 
ISmail et all 119971 : IBarger et alj I1998D . The first such 
galax ies were discovered via their submillimeter emis- 
sion (ISmail et al.lll997t IBarger et aflll998t iHughes et all 
1998), and are now co mmonly referred toas submillime- 
ter galaxies (SMGs). IChapman et alj (|2005f ) estimated 
that ~ 65% of such spectroscopically-confirmed bright 
SMGs (e.g., with Sgso/un ^ 5 mJy) at z ~ 2 — 3 have 
rest-frame UV colors similar to those of BXs and LBGs, 
yet are on average a factor of m 10 times more luminous. 
There is some uncertainty in the luminosities related 
both to the conversion between mid and IR luminosi- 
ties to total bolometric luminosities an d the fraction of 
the lu minosity that arises from an AGN (|Alexander et al.l 
I2005H . Taking the far-IR estimates of the SFRs of SMGs 
at face value then suggests that SMGs are examples of 
galaxies whose UV slopes typically under-predict their 
total attenuation and h ence total bolometric luminosi- 
ties (jReddv et al.ll2006b[ ). 

Measuring the frequency of such dusty galaxies among 
UV-faint sourc es requires th a t we estimate the former's 
space density. ICoppin et alj (|2006f ) determine a surface 
density of SMGs with ffssoit m > 5 mJy of 0.139 arc min~ 2 . 
The spectroscopic study of IChapman et all (|2005l ) found 
that 50% of bright SMGs lie at redshifts 1.9 < z < 2.7, 
implying a space density at these redshifts of 2.63 x 
10~ 5 ' Mpc~ 3 . These authors also found 30-50% of them 
have 25.5 < TZ < 28.0, corresponding to i uv < 0.34L* 
at the mean redshift of the BX sample (z — 2.30). 
According to our UV LF, the total number density of 
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galaxies over this same apparent magnitude range is 
3.28 x 1(T 2 Mpc -3 , implying that UV-faint SMGs with 
7Z > 25.5 constitute 0.02 — 0.04% of sources on the faint- 
end. Even in the most conservative case where we assume 
that all SMGs with Ssso^m > 5 mJy lie at 1.9 < z < 2.7 
and all have 1Z > 25.5, we fin d a fractional contribu tion 
of only 0.16%. The results of IChapman et all (|2005f) in- 
dicate that this SMG fraction would be even lower among 
z > 2.7 galaxies, although we note that their adoption 
of a radio-preselection may have biased the distribution 
of their sources to lower redshifts. We make note of the 
fact that the exact contribution will depend on the limit 
of what we consider to be "bright" SMGs, and extend- 
ing the limit to fainter submillimctcr fluxes will undoubt- 
edly include galaxies that are less attenuated, on average, 
and more likely to be recovered via their UV co lors (e.g., 
iReddv et all [20051 : lAdelberger fc Steidell l2000h . In any 
case, the current best estimates for SMGs that are ob- 
served routinely in the first generation of submm surveys 
imply that by number they make a very small contribu- 
ti on to the number de nsity of sub-L* galaxies. 

IReddv et alj (|2006bh demonstrate that the vast major- 
ity of luminous infrared galaxies (LIRGs) at z > 2 will 
have rest- frame UV colors that satisfy the BX/LBG cri- 
teria. While such criteria also pick up a non-negligible 
number of ultraluminous infrared galaxies (ULIRGs), 
the best method of accounting for these galaxies is via 
their infrared emission. The launch of Spitzer enabled 
observations that are sensitive to the warmer dust in 
high redshift starburst galaxies. Such galaxies are lu- 
minous in the infrared and appear to account for an in- 
creas i ng fraction o f galaxies at z > 1 (e.g.. iDev et all 
[20081 ICaputi et all [20071 iLeFloc'h et al.ll2005f) . Based 
on such studies, there have been a few estimates of 
the number density of ul traluminous infrared galaxies at 
z ~ 2 — 3. For instance, ICaputi et al.l (|2007l ) find a den- 
sity of (1.5 ±0.2) x 10~ 4 Mpc~ 3 for 24/zm-selected galax- 
ies with Lboi > 10 12 L (excluding AGN) in the GOODS 
fields. Similarly, 24/mi-bright galaxies (/24pm > 0.3 mJy) 
with red optical to mid-IR colors (R — [24] > 14) have 
a space density in the redshif t range 0.5 < z < 3.5 of 
(2.82 ± 0.05) x 10" 5 Mpc~ 3 ([Dev et al.1 [20081 ). almost 
all of which lie below L*, the characteristic unobscured 
UV lumin osity. Assuming con servatively that all of the 
ULIRGs of lCaputi etldl (|2007T ) are fainter than K = 25.5 
and as faint as 1Z ~ 28.0, we find a ULIR G fraction on the 
UV faint-end of 0.46%. In terms of the lDev et all (|2008l ) 
objects, assuming their space density does not evolve be- 
tween 0.5 < z < 3.5, this fraction is 0.086%. Hence, 
while such dusty, star-forming galaxies contribute signif- 
icantly to the total IR luminosity density, they must be 
vastly outnumbered by galaxies with fainter bolometric 
luminosities. This result is not surprising given the close- 
to-exponential drop-off in number counts of such infrared 
luminous galaxies according to the Schechter function, 
combined with the steep faint-end slope of the UV LF. 

Taken another way, if we make the supposition that 
a large fraction of galaxies on the faint-end of the UV 
LF are indeed very dusty star-forming ULIRGs, then 
by virtue of the sheer numbers of UV-faint galaxies, we 
would predict a number density of ULIRGs significantly 
in excess of the measured value. These calculations indi- 
cate that rapidly star-forming, dusty galaxies constitute 
a very small fraction of the total number density of star- 



forming galaxies on the faint-end of UV LF. Moreover, 
they support our premise that the E(B — V) distribu- 
tion is unlikely to be redder for UV-faint galaxies than 
for UV- bright ones (appendix). 



7.2. Galaxies with Large Stellar Masses 

Another population of galaxies at z > 2 charac- 
terized by their faint UV luminosities are those that 
have undergone their major episod e(s) of star forma- 
tion and are evolving quiescently (|Franx et alj [2003), 
commonly referred to as "Distant Red Galaxies," or 
DRGs . Such galaxies have low specific star fo rmation 
rates (|Papovich et al.ll2006t IReddv et al1l2006bD and are 
inferred to have low gas fractions ( Reddv et al.l l2006b) 
relative to UV-selected galaxies. The bulk of BX/LBGs 
have s tellar masses i n the r a nge 10 9 — 10 11 M^ (|Erb et all 
[20061 IReddv et al.l l2006at iShaolev et ail |2005j), while 
K < 20 DRGs have typical stel l ar ma sses of > 10 11 M 
(e.g., Ivan Dokkum et all 12004 I2006D . although there 
is some small overlap in the st ellar mass distribu- 
tion between BX/LBGs and DRGs (jShaplev et al.ll2005t 
Ivan Dokkum et al"1l2006lL part icularly at fainter iVband 
magnitude (IReddv et al1l2005l) . 

Ivan Dokkum et all ([2006) find that galaxies with stel- 
lar masses > 10 11 M Q are also typically faint in the op- 
tical, with > 2/3 fainter than 1Z = 25.5. Conservatively 
assuming that all such galaxies are fainter than 7Z = 25.5 
and as faint as 1Z ~ 28.0, and have an estimated space 
densi ty of (2.2 ± 0.6) x 10" 4 Mpc -3 (|van Dokkum et all 
l2006h . then we compute a fractional contribution to the 
faint-e nd of the UV LF in the sa me magnitude range of 
0.67%. Ivan Dokkum et alJ (|2006h noted that only 1/3 of 
these massive galaxies had the colors of BX/LBGs. How- 
ever, examination of their U n G7Z color distribution shows 
that a large fraction of the "missing" 2/3 have colors 
that hug the BX/LBG selection boundaries. Our incom- 
pleteness corrections will take into account objects that 
scatter into the BX/LBG samples because of stochastic 
effects like photometric errors, but conservatively assum- 
ing that 2/3 of massive galaxies are missed even after 
these corrections would imply a massive galaxy fraction 
among UV-faint sources of w 2%. 

These estimates imply that like dusty star-forming 
galaxies, those with large stellar masses (> 10 11 M Q ) 
comprise a very small fraction (< 2%) of all UV-faint 
galaxies. Hence, virtually all sub-L* galaxies have 
smaller stellar masses and are less dusty than the types 
of galaxies considered above. From a broader perspec- 
tive, several studies have shown that galaxies with large 
stellar masses ten d to cluster more s trongly than less 
massiv e galaxies (jQuadri et alj 120071 : lAdelberger et al.l 
l2005al) . This is consistent with the expectation that 
galaxies with large stellar masses formed stars earlier 
in more massive potential wells whic h are expected to 
be the most clustered. Furthermore, lAdelberge r et al.1 
(|2005d ) demonstrated that UV-bright galaxies cluster 
more strongly than UV-faint ones, at least at z > 2. 
Given the sheer number of UV-faint galaxies, these ob- 
servations suggest that galaxies on the faint-end of the 
UV LF are likely to be less clustered than their brighter 
counterparts, and hence associated with lower mass ha- 
los. 
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DISCUSSION: CONSTRAINTS ON THE STAR 
FORMATION RATE DENSITY 




Fig. 8. — Unobscured UV luminosity density, puv> per 0.5 
magnitude interval (dashed lines) and integrated (solid lines) at 
1.9 < z < 2.7 (cyan, blue) and 2.7 < z < 3.4 (magenta, red), 
respectively. Dotted lines indicate M* at 2 ~ 2 a nd z ~ 3. 
The e quivalent star formation rate density assuming the K ennicu"t9 
( 1998) relation and a Kroupa IMF is shown on the right-hand axis. 



As is customary, the iKennicuttj (|1998f ) relation is used 
to convert UV luminosity to star formation rate (SFR), 
adopting a Kroupa IMF from 0.1 to 100 M Q (Figure [SJ 
Table H Table (3). This results in factor of - 1.7 de- 
crease in SFR for a given luminosity owing to the larger 
fractional cont ribution of hig h-mass stars to the Kroupa 
relative to the lSalpeterl (1955) IMF. For consistency with 
previous investigations, the luminosity density is calcu- for extinction using a method similar to that presented 
lated to a limiting luminosity of 0.04L* =3 unless stated in R08. Briefly, a large number of galaxies are simu- 



UV colors — as motivated by S pitzer mid-IR obse rva- 
tions of UV-selected g alaxie s (|Reddv et al.l l2006bl ) - 
via the iCalzetti et al.l (|2000T ) relation. This has been 
shown t o be valid for moderately lu minous galaxies (i.e., 
LIRGs; IReddv et al.l I2005L I2006bf) . We will consider 
shortly the contribution from high redshift galaxies that 
do not follow the local starburst atten uation relations 
(jMeurer et al.lll999fc ICalzetti et al.ll2000fl . 

Defining N{E{B - V),L) = N(E(B - V)) will of 
course leave the relative contribution of UV- faint galaxies 
to pjjv unchanged from their contribution to the unob- 
scured UV luminosity density. The bolometric luminos- 
ity density is calculated under the more realistic case of a 
declining average reddening with unobscured luminosity 
(appendix), with the results tabulated in Tables [4] and OH 
In this case, we find that galaxies fainter than L* — de- 
fined as the characteristic unobscured UV luminosity — 
contribute 0.62 and 0.78 at z ~ 2 and 3, respectively, to 
Puy integrated to 0.04L* =3 . These fractions are likely to 
be lower limits since there is a non-negligible number of 
very dusty and bolometrically luminous UV-faint galax- 
ies at these redshifts (§[7|). Below, we revisit our estimate 
of the bolometric luminosity density after incorporating 
the effect of the most luminous galaxies at z ~ 2 — 3. 

8.2. Bolometric Luminosity Functions 

Although so far we have defined L* in terms of the 
knee of the UV LF uncorrected for extinction, we can 
also examine the fractional contributions as a function 
of luminosity to the bolometric luminosity density. This 
is accomplished by reconstructing the UV LF corrected 



otherwise. The differential and cumulative unobscured 
UV luminosity densities to 0.04L* =3 , p\jv(> 0.04L* =3 ), 
are « 6% and 42% larger at z ~ 2 and 3, respectively, 
than reported by R08. This difference is attributable to 
the steeper faint-end slope and slightly brighter L* de- 
rived in this study. Below, we consider the effects of a 
luminosity-dependent dust correction, the bolometric lu- 
minosity functions at z ~ 2 — 3, and implications for the 
star formation history. 

8.1. Luminosity- Dependent Dust Corrections 

As a consequence of the steep faint-end slopes at 
z ~ 2 — 3, ss 93% of the unobscured UV luminosity 
density (integrated to zero luminosity) is contributed by 
galaxies fainter than L* (Figure [5]) . The abundance of 
UV-faint galaxies and their cumulative luminosity makes 
them ideal candidates for the sources responsible for most 
of the ionizing flux at z > 3. However, the luminosity 
dependence of reddening implies that their contribution 
to the bolometric luminosity is likely to be diminished 
compared to their contribution to the unobscured lumi- 
nosity density. The bolometric luminosity density can be 
expressed simply as 



bol 



L<t>(L)10l° Ak 'W A WdL, 



(3) 



where A(L) is the reddening, parameterized by E(B — V), 
as a function of luminosity and fc'(A) is the s t arburs t 
attenuation relation defined in Calzet tTet al.l ((2000) . 
For this calculation, we have assumed that the bolo- 
metric luminosity can be recovered from the rest-frame 



lated with magnitudes and E(B — V) drawn randomly 
from the LF and luminosity-depe ndent E(B — V) distri- 
bution. The Caiz etti et al.l (|2000f ) relation is used to re- 
cover the bolometric luminosities, which are then binned 
to produce a luminosity function (Figure [§])■ We allow 
for the LF to vary within the errors and add a 0.3 dex 
scatter to the dust correction implied by E(B — V), re- 
flectin g the approximate scatter in both the local rela- 
tions JMeurer et al.lll999tlCalzetti et al.ll2000f ) and those 
found at z ~ 2 ( Reddv et al. 2006b). This scatter results 
in a 5% random error in the faint-end of the bolometric 
LF, significantly smaller than the systematic errors that 
result from assuming different relations between dusti- 
ness and UV luminosity (Figure [9]). 

Note that we use only the E(B — V) distribution found 
for UV-selected galaxies to reconstruct the bolometric lu- 
minosity functions. The range in attenuation factors ob- 
tained for such galaxies will be smaller than the intrinsic 
range of reddening among all galaxies at z ~ 2 — 3. One 
obvious reason for this bias is the incompleteness for ob- 
jects that never scatter into our sample because of their 
red colors. Another reason is that even if such red, dusty 
galaxies do satisfy the LBG color criteria, their bolomet- 
ric luminosities may be underestimated severely based 
on the UV colors alone. Hence, the method for recover- 
ing bolometric LFs based on the E(B — V) distribution 
of galaxies that scatter into the BX/LBG windows will 
underestimate the contribution of galaxies to the bright- 
end of the bolometric luminosity function. Because of 
this, the contribution of these dusty galaxies is based on 
results published elsewhere. Specifically, we adopt the 
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TABLE 4 

Total UV Luminosity Densities at 1.9 < z < 3.4 



Redshift Range 


Unobscured 3, 
(ergs s" 1 Hz -1 Mpc -3 ) 


Dust-Corrected b 
(ergs s" 1 Hz -1 Mpc -3 ) 


1.9 < z < 2.7 
2.7 < 2 < 3.4 


(3.89 ± 0.24) x 10 26 
(3.28 ±0.24) x 10 20 


(1.36 ±0.30) x 10 27 
(8.74 ±2.55) x 10 26 



a Uncorrected for extinction, integrated to 0.04 L*_ 3 . 
b Corrected for luminosity-dependent extinction, including both ob- 
scured and unobscured UV luminosity, integrated to 0.04 L*_ 3 . 



TABLE 5 

SFRD Estimates and Dust Correction Factors 





2 r 


- 2 




- 3 




Llim = 0.04L* a 


L lim = 


L lim = 0.04L* a 


Llim = 


(1) UV SFRD uncor b 


0.032 ± 0.002 


0.064 ± 0.003 


0.027 ±0.002 


0.055 ±0.003 


(2) UV SFRD cor (LDR) bc 


0.112 ±0.025 


0.122 ±0.027 


0.072 ± 0.021 


0.080 ± 0.023 


(3) UV Dust Correction (LDR) C 


3.50 ±0.78 


1.91 ±0.42 


2.67 ± 0.78 


1.45 ±0.42 


(4) UV SFRD cor (CR) bd 


0.144 ±0.009 


0.288 ±0.014 


0.122 ±0.009 


0.248 ± 0.014 


(5) UV Dust Correction (CR) d 


4.5 


4.5 


4.5 


4.5 


(6) Total SFRD cor bc 


0.142 ±0.036 


0.152 ±0.038 


0.102 ±0.032 


0.110 ±0.034 


(7) Total Dust Correction 1 


4.44 ± 1.13 


2.38 ±0.59 


3.78 ± 1.19 


2.00 ±0.62 



a Integrated to include all galaxies with unobscured UV luminosities brighter than 0.04L*_ 3 , or M(1700A) ^ —17.48. 

b In Mq yr _1 assuming a Kroupa IMF. 

c Invokes luminosity-dependent reddening (LDR). 



d Invokes luminosity-invariant reddening (CR). 
e Sum of LDR-corr ected star formation rate density from row(2) and the contribution of -Lboi > 10 12 L0 galaxies from 
ICaputi et all pOOl ). 

f Dust correction required to recover total star formation rate density in row (6) from the unobscured star formation rate 
density in row (1). 

track the bright-end of the IR LF. The results are sum- 
marized in Figure [9] assuming no evolution in the bright- 
end of the LF (Lboi ^ 10 12 L Q ) in the redshift range 
1.9 < z < 3.4. It is important to keep in mind that 
vL y at 1700 A scales with SFR in a different way than 
the infrared luminosity (Lir = L(8 — 1000//)). Hence, 
the bolometric luminosity — the sum of the UV and IR 
luminosities as defined in this paper — will scale in a 
non- linear way with SFR. In the present context, star 
formation rate densities are computed separately for (1) 
galaxies where Lboi is determined from the UV-corrected 
values and (2) galaxies with Lboi ^ 10 12 L Q where the 
bolometric lumino sity is determined p urely from the in- 
frared luminosity (|Caputi et al.|[2007l ). The star forma- 
tion rate densities from the two contributions are then 
added to estimate the total. With the appropriate scal- 
ings, this calculation implies that ~ 70 — 80% of the 
bolometric luminosity density arises from galaxies with 
Lho\ 10 12 L Q , consistent with findings of R08. 

Taken together, these findings can be summarized as 
follows. Including ULIRGs — those galaxies whose UV 
slopes te nd to under-predict their bolometric luminosi - 
ties (e.g.. lPapovich et al.ll2006ll2007l : lReddv et al.ll2006bl ) 
— does not change the fact that a large portion of the 
bolometric luminosity density arises from faint galax- 
ies, either those that are fainter than the characteris- 
tic unobscured UV luminosity or those that are fainter 
than the characteristic bolometric luminosity. Placing 
these results in a wider context will require more precise 
estimates of the bright-end of the bolometric luminos- 
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Fig. 9. — Bolometric luminosity functions at 2 ~ 2 (blue) and 
2 ~ 3 (red), computed by combining the measurement of the UV 
luminosity function with a luminosity-dependent E(B — V) dis- 
tribution (see text). The upper limits of the shaded regions indi- 
cate the LF derived assuming a constant E(B — V) distribution. 
The lower limits indicate the LF derived assuming that all galaxies 
with apparent magnitude fainter than 1Z = 25.5 have zero redden- 
ing. These limits encompass the range of likely LFs and give an 
indication as to the systematic uncertainty in the bolometric LF. 
The solid lines denote the bolometric LF obtained using our model 
of the luminosity-dependent E(B — V) distribution that gradu- 
ally falls to zero reddening for the faint est galaxies. At 2 ~ 2, 
the higher luminosity points (circles) from Caputi et al. ( 20071) are 
shown, along with their Schechter extrapolation to fainter lumi- 
nosities (long-dashed line). 

value of the bright-end of the infrared LF at z ~ 2 (af- 
ter ex clusion of bright AGN) presented by ICaputi et al.1 
(|200l since the bright-end of the bolometric LF should 
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ity function that (1) take into account the luminosity- 
dependent conversion between mid-IR luminosity (upon 
which most estimates are based) and the total infrared 
luminosity and (2) the potential contamination from 
AGN that are prevalent among galaxies with such high 
IR luminosities. Nonetheless, combining the most re- 
cent estimate of the bright-end of t he bolometric lumi- 
nosity function (jCaputi et al.ll2007l ) with our results at 
the faint-end points to a luminosity density that is dom- 
inated by bolometrically faint to moderately luminous 
galaxies. The implications for a luminosity-dependent 
reddening distribution on the average dust correction fac- 
tors applied to high redshift samples and the evolution 
of the star formation rate density are discussed in the 
following sections. 

8.3. Average Dust Correction Factors 

A luminosity dependent dust correction and the large 
ratio of UV-faint to UV-bright galaxies implies an aver- 
age UV dust correction that is sensitive to the limit of 
integration used to compute the luminosity density. It 
seems prudent to consider such a systematic effect given 
that estimates of the star formation rate density imply 
stellar mass densities in ex cess of what are actually mea- 
sured (Wilkins et al.l 120081 ). This effect is mentioned in 
R08; here, we proceed to quantify the average dust cor- 
rection factors relevant for luminosity densities computed 
to different limits based on our new determination of the 
faint-end slope. 

The calculated dust corrections and star formation rate 
densities are listed in Table [5l We have assumed a contri- 
bution of Lboi > 10 12 L Q galaxies to the s tar formation 
rate d ensity at z ~ 2 as computed from ICaputi et al.1 
([2001 . We also assume this same contribution at z ~ 3, 
though it has not been measured directly at these higher 
redshifts, in order to place conservative estimates on the 
effect of a luminosity-dependent dust correction on the 
average dust correction factors required to convert UV 
luminosity densities to star formation rate densities. 

The luminosity-dependent reddening model implies 
dust corrections of a factor of 3.5 and 2.7 at z ~ 2 and 3, 
respectively, integrated to 0.04L* (TableEJ), which are up 
to a factor of two smaller than the typical 4.5 — 5.0 dust 
corrections found for 1Z < 25.5 galaxies (jSteidel et al.1 
119991: iReddv fc Steidelll2004h . Aside from differences in 
the luminosity range probed, this difference in average 
extinction is mitigated somewhat by the fact that a sig- 
nificant fraction (~ 0.2 — 0.3) of the bolometric lumi- 
nosity density arises from ULIRGs, where the usual dust 
conversions do not apply (see discussion above). The 
expectation is that the lower dust-corrected luminosity 
densities inferred in the luminosity-dependent reddening 
case are compensated by the inclusion of galaxies where 
E{B — V) tends to under-predict the reddening. The to- 
tal dust corrections required to recover the bolometric lu- 
minosity density, including that contributed by ULIRGs, 
are 4.4 and 3.8, somewhat larger than the values quoted 
above. While the differences between these dust correc- 
tions may seem small at z ~ 2 — 3, they do result in 
up to a factor of two in systematic scatter in star for- 
mation rate density measurements, comparable to the 
dispersion in the local calibrations between luminosity 
and star formation rate, and so should be taken into ac- 
count. The dependency of the average dust correction 



on the integration limit will be even greater for steeper 
faint-end slopes given the larger fractional contribution 
of less-reddened faint galaxies to the luminosity density. 
The average dust correction factors stated above are rel- 
evant when integrating the UV luminosity function to 
0.04L* =3 . Integrating to zero luminosity alters the cor- 
rections to be a factor of 2.4 and 2.0 at z ~ 2 and 3, 
respectively. To reiterate, these extinction corrections 
account for not only the dust-obscuration among moder- 
ately luminous galaxies prone to UV-selection, but also 
for those ultraluminous galaxies that may either escape 
UV-selection or simply have rest-UV slopes that under- 
predict their bolometric luminosities. These effects un- 
derscore the various subtleties that can affect extinction 
corrections for UV-selected samples. 

Of course, these dust corrections are equally impor- 
tant at higher redshifts z > 3 where the only constraints 
on the star formation rate density come from UV ob- 
servations. Evidence suggests that UV-selected galaxies 
become bluer a t redshifts z > 3, relative to galaxies at 
lower redshifts (|Yan fc Windhorst! 120041 : iBouwens et al.1 
I2007I) . 6 This trend may be attributable to two effects. 
First, as noted above, L* evolves strongly as a func- 
tion of redshift at z > 3, such that the average UV lu- 
minosity of galaxies decreases with increasing redshift. 
If sub-L* galaxies have lower dust reddening than UV- 
bright ones, the trend in UV color may be interpreted 
as a decrease in dust reddening. A consequence of the 
luminosity-dependent reddening model is that, when ex- 
amined over a large dynamic range of luminosity, the 
unobscured UV luminosity must track the bolometric lu- 
minosity, or SFR. This leads to the second effect which 
is tied to the observation that high redshift galaxies are 
less attenuated than lower redshift galaxies of the same 
bolometric luminosity, resulting in a trend of decreas- 
ing ex tinction per unit SF R proceeding to higher red- 
shifts (|Reddv et al.l I2006E R08). Stated anoth er way, 
extrapolating the results of IReddv et al.1 (|2006b( ) to red- 
shifts z > 4 implies that higher redshift galaxies are less 
opaque at UV wavelengths than lower redshift galaxies 
of the same bolometric luminosity. Hence, the combined 
effect of a lower L* and lower average dust attenuation to 
a given bolometric lu minosity implies lower dust correc- 
tions at higher z (R08. IBouwens et al.ll2007t l. The impli- 
cations for this evolution in the average dust correction 
are discussed in the next section. 

There have been several explanations put forth to ex- 
plain the discrepancy between the integrated star for- 
mation history and stellar mass density measurements, 
includ i ng m i ssing stellar mass, a n evolution of the IMF 
(| Pavel 120081 : Ivan Dokkuml [2008) , or more generally an 
evolving conversion between luminosity and star forma- 
tion rate. In light of these effects, we consider the po- 
tential impact of an evolving dust correction on the star 
formation history, as described below. 

9. DISCUSSION: STAR FORMATION HISTORY AND 
BUILDUP OF STELLAR MASS 

e It is generally accepted that the highest redshift star-forming 
populations exhibit bluer UV colors on average than similarly 
selected galaxies at lower redshifts, although the result has not 
been verified indepen dently with longer wavelength observations. 
IBouwens et al.1 1120071 ) present a discussion of why there may not 
be a strong bias against dusty galaxies with the higher redshift 
dropout criteria. 
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We have already touched upon a few of the implica- 
tions of a steep faint-end slope on the star formation rate 
density. In particular, we noted that a very large frac- 
tion of the unobscured UV luminosity density (> 90%) 
arises from galaxies fainter than the characteristic unob- 
scured UV luminosity. Similarly, our results suggest that 
even assuming a lower reddening among UV-faint galax- 
ies relative to UV-bright ones implies a bolometric, or 
dust-corrected, luminosity density dominated by galax- 
ies fainter than the characteristic bolometric luminosity. 
There are several important consequences of these results 
that we discuss in the next few sections. 

9.1. Contributions at z ~ 2 — 3 to the Global Stellar 
Mass Density 

A significant fraction of the stellar mass density that 
formed between z — 1.9 and z — 3.4 (the redshift limits 
of our analysis) — corresponding roughly to the epoch 
when galaxies were forming most of their stars (§ [JJ - 
occurs in galaxies with Lboi ^ 10 12 L Q . Using a linear 
interpolation of the contributions of galaxies with dif- 
ferent luminosities to the bolometric luminosity density 
between z ~ 3 and z ~ 2 and multiplying by the time 
between z = 3.4 and z — 1.9, yields a total SMD of 
0,(1.9 < z < 3.4) = 0.0014 ± 0.0003 in units of the 
critical density. This value is alre ady 0.57 times th at of 
the present-day value reported in ICole et alJ (|2001h . As 
mentioned above, there is still a fair amount of uncer- 
tainty regarding the bright-end of the bolometric lumi- 
nosity function. Irrespective of the number density of 
bolometrically-luminous galaxies, our calculations sug- 
gest that 43 ± 9% of the present-day stellar mass density 
was formed in galaxies with 6 x 10 s < Lboi < 10 12 L Q be- 
tween redshifts z = 3.4 and 1.9. 7 While much attention 
recently has been focused on the stellar mass buildup 
associated with luminous galaxies at high redshift, it is 
clear that fainter galaxies, those that are routinely picked 
up in UV surveys but may be missing from rest-optical 
and far-IR ones, also play an important role. It suggests 
that much of the stellar mass assembly at a time when 
galaxies where forming most of their stars occurred in 
the typical and more numerous galaxies that populate 
these redshifts. 

9.2. Evolution of the Star Formation Rate Density 

A luminosity-dependent dust correction not only has 
consequences for the total star formation rate density 
measured at a given redshift, but the strong evolution of 
L* suggests that it will induce a systematic effect with 
redshift. Figure [TOl summarizes the star formation rate 
densities inferred with the luminosity-de pendent redden- 
ing m odel from this study and those of iBouwens et all 
(2007) at z > 3.8, compared with the star formation 
history assuming a constant dust-correction of 4.5 (red 
li ne), along with lower redshift determinations compiled 
in Hopkins] (|2004h . For consistency with the latter study, 
we have integrated to zero luminosity (see § 19 .41 for a dis- 
cussion of the systematic effects associated with the lim- 
its of integration). The luminosity-dependent reddening- 
corrected star formation history points to a factor of 8 — 9 

7 The limit of 6 X 10 8 Lq is adopted for consistency with R08. 
The bolometric LF exhibits a slope that is somewhat shallower than 
the UV LF, so changing the limit of integration to zero bolometric 
luminosity will add roughly 10% to the luminosity density. 



increase in the star formation rate density between z ~ 6 
and z ~ 2 (e.g., see also R08. IBouwens et al"1l2007| ). sig- 
nificantly steeper than the factor of 4 that we would have 
inferred in the case of a luminosity-invariant (constant) 
dust correction. As discussed in § 16.11 the evolution in 
the unobscured star formation rate density is connected 
to the increase in number, and hence luminosity, density 
of galaxies brighter than L*. Figure [TOl demonstrates 
that the comparable evolution in the bolometric star for- 
mation rate density is driven both by an increase in the 
number density of bright galaxies and an evolving dust 
correction. 

The elevated star formation rate densities predicted 
in the luminosity-invariant reddening model result in a 
stellar mass growth of O = 0.0026 ± 0.0007 between 
2.3 < z < 5.9, compared to O = 0.0013 ± 0.0003 
for the luminosity-dependent reddening model over the 
same redshifts, a factor of two difference between the two 
dust correction scenarios. More importantly, a constant 
dust correction model predicts stellar mass buildup be- 
tween 2.3 < z < 5.9 that exceeds the local measurement. 
Notwithstanding the noted disagreement between stel- 
lar mass density measurements and the integral of the 
star formation history (Figure [TU] and discussion below) , 
these results suggest that we may be able to rule out the 
elevated star formation rates predicted by constant dust 
correction models, although we caution that the differ- 
ences in total stellar mass accumulated by z ~ 2.3 based 
on a constant versus declining star formation history are 
small compared to the uncertainties in the star forma- 
tion rate and stellar mass density measurements. In the 
next section, we discuss this disparity between the in- 
tegrated star formation history and stellar mass density 
measurements and possible resolutions. 

9.3. Reconciling the Star Formation History with 
Stellar Mass Density Measurements: 
Luminosity-Dependent Dust Corrections and 
Missing Stellar Mass 

Figure [10] shows the star formation history inferred 
by differentiating measurements of the stellar mass den- 
sity (integrated to zero) as a function of redshift (purple 
hashed region). The results imply that there is a max- 
imum disparity of w 0.5 dex in this inference and ac- 
tual observations of the star formation rate density at 
z ~ 2 — 3. It is of general interest to determine whether 
this discord is due to some lack of understanding of the 
fundamental physical processes that govern star forma- 
tion and/or to the mundane nature of the uncertainties 
that seemingly plague SFR and stellar mass estimates, 
including sample incompleteness and the limits to which 
one integrates to obtain the star formation rate and stel- 
lar mass densities. 

In light of the steep faint-end slopes of the UV LF 
advocated at z > 2, it is worthwhile to consider the pos- 
sibility that the stellar mass density measurements at 
these redshifts are too low, primarily because they do 
not account for low mass galaxies that may escape stel- 
lar mass selected samples but, even with their low stellar 
masses, are sufficiently numerous to add appreciably to 
the total budget of stellar mass. The comparison drawn 
in Figure [10] implicitly assumes that all the galaxies con- 
tributing to the estimate of the star formation rate den- 
sity are in some way also represented in the estimate 
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Fig. 10. — Cosmic star formation history, including the 
luminosity-dependent dust-corrected determinatio ns at z ~ 2 — 3 
from this analysis (large pentagons) and those of Bouwcns ct al. 
j2007h at z > 3 .8 ( open circles at high z), and the compilation from 
iHopkinsi mm ( open squares) at low z. Note that our estimates 
include the directly-measured contribution to the star formation 
rate density from ultraluminous infrared galaxies and assume that 
this contribution is non- evolving between z ~ 3 and z ~ 2. Also, 
for consisten cy with th e Hopkins (2004) compilation, our points 
and those of Bouwcns ct al. I 1 120071) are computed by integrating 
the UV LF to zero luminosity. The solid red line shows the star 
formation history assuming a constant dust correction of 4.5 to 
the unobscured UV luminosity densities at z ~ 2 — 6. The short- 
dashed line shows the fit to the star formation history including this 
constant dust correction model. The solid blue line indicates the 
best-fit star formation history assuming a luminosity-dependent 
dust correction to the z > 2 measurements. The solid hatched pur- 
ple region indicates the ±1 a star form ation history inferred from 
the evolution of the stellar mass density (jWilkins et al.l 2003) , with 
an extrapolation at z > 4.5 based on st e llar mass density measure - 
ments at z ~ 5 - 6 fromlMcLure et alj <200lj); lEvles et al.l ||2007T) ; 
IStark et all (I2007T ) : IVerma et al.l 1120071) ; mm et al.l ( dashed 

purple region). As discussed in the text, much of the discrepancy 
between the stellar mass density measurements and the integral of 
the star formation history may be due to incompleteness of low 
mass objects in the stellar mass estimates. A Kroupa IMF is as- 
sumed throughout. 

of the stellar mass density. In practice, the problem is 
that unlike SFR- limited samples, mass-selected samples 
at high redshift do not probe far enough down the stel- 
lar mass function due to the significant amount of time 
required to assemble the requisite near-IR data. Hence, 
such studies may underestimate the low-mass slope of 
the stellar mass function. 

9.3.1. Stellar Mass Density in UV-Bright (TZ < 25.5) 
Galaxies 

The slope of the stellar mass function at z ~ 2 — 3 
is not well-constrained. However, if we are able to es- 
timate the average stellar mass of LBGs, then knowing 
their number density from the UV LF will enable us to 
estimate their contribution to the stellar mass density. 
We compiled stellar mass estimates for BXs and LBGs 
in the GOODS-N an d Q1700 fields (|Reddv et al J 120061 : 
Shapl ev et al.l 12001 ). Briefly, stellar masses are com- 
puted for spectroscopically-confi rmed BXs and LBGs by 
fitting iBruzual fc Chari ot (2003) model templates to the 
observed U n GlZ+ J K S +IRAC photometry, and allowing 
the star formation history r and E(B — V) to vary freely. 
The star formation rate and stellar mass are determined 
by the normalizati on of the model SEP to the broad- 
band photometry jShaplev et al.ll2005t lErb et "ail 120061 : 
IReddv et al.ll2006aft . 

Excluding spectroscopically identified AGN, the distri- 
butions of stellar mass for 208 BXs and 42 LBGs from 



the 2 aforementioned fields are shown in Figure fTTl Since 
most moderately star-forming galaxies escape BX/LBG 
selection due to stochastic effects (e.g., photometric scat- 
ter; §|3l), we adopt the reasonable premise that galaxies 
to TZ = 25.5 that do not satisfy the BX/LBG criteria 
have a similar distribution in stellar masses to those that 
do. Note that, as discussed previously, this may not be 
the c ase for the most massive galaxies at these redshifts 
(e.g.. Ivan Dokkum et al.ll2006h . However, in the present 
context we are interested primarily in the contribution 
of typical star-forming galaxies — which outnumber by 
far the most massive galaxies at these redshifts (§ EJ - 
to the stellar mass density. Further, we show in the ap- 
pendix that adopting a young stellar population does not 
affect appreciably the incompleteness corrections to the 
BX/LBG samples. Consequently, it is unlikely that large 
numbers of young galaxies with low stellar masses are 
scattered out the sample relative to the frequency with 
which more massive galaxies are scattered out of the sam- 
ple, particularly among UV-bright galaxies. Therefore, 
we make the reasonable assumption that the stellar mass 
distributions for UV-bright BX/LBGs are representative 
of UV-bright star-forming galaxies in general. Figure [TT] 
also shows the distribution in absolute UV magnitude 
for galaxies with stellar mass estimates, spanning the 
full range of magnitudes re presented in the spec troscopic 
sample. As first shown bv lShapley et al.l (|2005h . we find 
no significant correlation between unobscured UV mag- 
nitude and stellar mass, perhaps not surprising since the 
two quantities are related only peripherally. We will re- 
visit this issue below. 

Based on these distributions, let us proceed to esti- 
mate the stellar mass density contributed by UV-bright 
galaxies. To do this, we generated many random re- 
alizations of the UV LF as allowed by the errors, and 
drew random absolute magnitudes from each of these re- 
alizations. We assigned a stellar mass drawn randomly 
from the observed distribution ( Figure [TTj) , which is then 
pert urbed by 0.3 dex to account for random uncertain- 
ties (jShaolev et al.l 120051) . The masses are then binned 
to produce a rough proxy for the stellar mass func- 
tion. The resulting Gaussian distributions for TZ < 25.5 
galaxies at z ~ 2 and 3 - corresponding to galaxies 
with Af A B(1700l) < -19.53 and -20.05 and z = 2.3 
and z — 3.05, respectively — are shown in Figure [T21 
Integrating these distributions for those galaxies with 
M* < 10 11 M Q yields: 

n*(M AB (1700A) < -19.53; < 10 11 M Q ; z ~ 2.3) = (4) 

(3.72 ±0.28) x 10~ 4 

Q*(M A b(1700A) < -20.05; < 10 11 M ; z ~ 3.05) = 

(1.53 ±0.15) x 10~ 4 (5) 

in units of the critical density (Table [6]). These estimates 
are meant to reflect the stellar mass densities contributed 
by UV-bright star-forming galaxies. Note that the stellar 
mass densities computed here differ from those derived in 
§ 19. 1( the latter are based on integrating the star forma- 
tion rate density, whereas the former are based on masses 
determined from broadband fitting of galaxy SEDs, and 
so are subject to somewhat different systematics. The 
important result of this section is that even without cor- 
rections for (1) the most massive and dusty galaxies at 
these redshifts for which the BX/LBG criteria are incom- 
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Fig. 11. — (Left): Distribution of stellar mass for spectroscopically-confirmed BX/LBG galaxies, excluding AGN, with mean values 
indicated in the panel and by the vertical lines. (Right): Stellar mass as a function of unobscured absolute UV magnitude for BX galaxies 
(filled squares) and LBGs (open circles). 



TABLE 6 

Stellar Mass Density Budget at 1.9 < z < 3.4 





1.9 < z< 2.7 


2.7 < z < 3.4 


n t (Tl < 25.5; < 10 11 M Q ) a 


3.72 ±0.28 


1.53 ±0.15 


n*(72. > 25.5; < 10 11 M ) b 


2.06 ± 0.26 


1.86 ±0.21 


n,(> io 11 m q ) c 


1.64 ±0.45 




Q„(Total) d 


7.42 ± 0.59 


> 3.39 ±0.26 



a Stellar mass density, in units of the critical density xl0 — 4 , 
in galaxies with TZ < 25.5 and stellar masses < 10 11 Mq 
assuming a Kroupa IMF. 

b Same as (a), but includes the contribution inferred for 
1Z > 25.5 galaxies based on the correlation between SFR 
and stellar mass for UV-selected galaxies (see text). 
c Stellar mass density in gal axies with stellar masses > 
10 11 Mq , based on the data of Ivan Dokkum et~aTI {2006). 
d Total stellar mass density, computed by adding the num- 
bers from the first three rows, including the contribution of 
UV-bright and faint galaxies, as well as those with stellar 
masses > 10 11 Mq. 

plete (§ [7]) and (2) UV-faint galaxies with 1Z > 25.5, we 
already find a stellar mass density at z ~ 2 comparable 
to estimates from rest-frame op t ically-selected samples 
("Figure f!2l iFontana et all 120031: [Dickinson etai] l2003t 
iRudnick et alj l2003t iDrorv et all l20"05h . Of course, di- 
rect comparisons between our measurements and those 
from optically-selected samples are fraught with signifi- 
cant biases, both random (e.g., field-to-field variations in 
the optical samples) and systematic (adopted rest-frame 
optical limits, underestimates of stellar mass by assuming 
a single component SF model, or more generally system- 
atics in the assumed M/L ratio and differences in stellar 
population models). Some of the random uncertainties 
are constrained by taking values from different surveys 
conducted in spatially disjoint fields, and at face value, 
the results above suggest that typical star-forming galax- 
ies already contain an amount of stellar mass comparable 
to that detected in rest-frame optically-selected surveys. 

9.3.2. Massive Galaxies 

From the survey results of Ivan Dokkum et all (|2004f ). 
after converting to a common IMF, the mass density con- 
tributed by galaxies with stellar masses > 10 11 M at 



2.0 < z < 3.0 is n*(> lO n M ) = (1.64 ±0.45) x 10" 4 
(Table [6]), where the uncertainty does not include po- 
tentially la rge systematic error s in photometric redshifts 
(e.g., R08. IShaplev et al1l2005D . For this calculation, we 
have assumed that the mass density does not evolve over 
redshifts 2.0 < z < 3.0, although it most likely does, 
and have assumed the aforementioned value is valid at 
z ~ 2.3. Adding this to the contribution from UV- 
bright galaxies yields a mass density of f2*(^ ~ 2.3) = 
(5.36±0.53) x 10 -4 . Do UV-faint galaxies contain enough 
stellar mass to add appreciably to this number? We ex- 
plore this question in the next section. 

9.3.3. Stellar Mass in UV-faint Galaxies 

How might the stellar mass distributio n be expected 
to cha nge for UV-faint s ub-L* galaxies? iShaplev et al.l 
(120051 ) and lReddv et alj (|2006bD highlight the biases in- 
herent in photometric redshift estimates for star-forming 
galaxies at z ~ 2, perhaps even more so for UV-faint 
galaxies. We are at present targeting UV-faint galaxies 
with deep spectroscopy to remedy this situation. Such 
spectroscopy, combined with deep multi- wavelength data 
in several of our survey fields, should allow us to con- 
strain the stellar populations and masses of sub-L* to as 
much confidence as one can obtain with such an analy- 
sis. A full SED analysis of such galaxies is beyond the 
scope of this paper, yet we can make some progress in 
determining the stellar mass content of sub-L* galaxies 
based on observations of UV-bright galaxies. 

To do this, we exploited the log- linear relation be- 
tween SFR and stellar mass fou nd at z ~ 2 . 3 from 
an analysis of deep HDF data bv ISawicki et al.l (|2007h : 
log(M stars /M ) = 9.0 + 0.86 log[SFR/(M e yr-% They 
present evidence that this correlation remains valid 
for galaxies with SFRs of w 1 M yr -1 , correspond- 
ing to unobscured UV magnitudes of Mab(1700A) ~ 
— 18.0. Several other correlations between SFR and stel- 
lar mass have been publishe d, includin g most recently by 
iReddv et all (|2006bD and lDaddi et~all (|2007f ). Adopting 
these latter relations results in a slightly larger contri- 
bution of stellar mass from UV-faint galaxies. There- 
fore, as a conservative estimate, we have adopted the 
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Fig. 12. — (Left): Stellar mass functions for 1Z < 25.5 star-forming galaxies at 1.9 < z < 2.7 (filled circles) and 2.7 < 2 < 3.4 (open 
squares), based on combining the number density computed from the UV LF and the stellar mass distribution measured for BXs and LBGs 
(Figure [TT} . The dotted lines indicate the inferred contribution from galaxies fainter than 1Z = 25.5, based on the trend between SFR and 
stellar mass for UV-selected galaxies (see text). The dashed lines indicat e the total c o ntribu tion from both UV- bright and faint galaxies. 
For comparison, the GOODS and FDF results at 2.25 < z < 3.00 from IDrorv et al.l | |2u05|) a re denoted by t he ope n triangles. ( Right) : 
Stellar mass d e nsity measurements at 1.9 < z < 2.7 (open triangles) f rom the following sou rces: IRudnick et~aTl p003 1 : IDrorv et alj <2005): 
IPozzetti et al.l H2007I 1: IFontana et al.l (120031 1: IDickinson et al.l [|2003I 1: IFontana et al.l I|2006I 1. All of these studies constrain the SMD over 
areas that are significantly smaller than the almost 1 square degree probed in this study, and most rely on photometric redshifts. Our 
estimates at z ~ 2.3 are shown by the large circles: dotted shows the estimate for UV-bright (7Z < 25.5) galaxies with M 4 < 10 11 Mq; 
dashed shows the estimate including UV-faint galaxies with M, < 10 11 Mq to a faint limit of Mab(1700A) = —18.0; and solid denotes the 
total co ntri bution including massive (NL > 10 11 M Q ) galaxies. The stellar mass density inferred by integrating the star formation history 
(Figure llOt to a limit of 0.083L* (same as that used to compute the stellar mass density) is denoted by the solid line. For comparison, the 
dashed line shows the result when integrating the star formation history to zero luminosity. 



ISawicki et al.l (|2007l ) correlation in the subsequent dis- 
cussion. 8 Although we do not observe a significant cor- 
relation between unobscured UV luminosity and stellar 
mass for UV-bright (7Z < 25.5) galaxies, the log-linear 
behavior of SFR with stellar mass implies that such a 
correlation must exist when examined over a large dy- 
namic range in unobscured UV luminosity. In particular, 
since UV-faint galaxies are likely to have lower reddening 
than their brighter counterparts (§0, the UV luminos- 
ity for these galaxies is expected to track the bolomctric 
luminosity gi yen the tight relation between SFR and red- 
dening (e.g.. iReddv et a l. 2006a). This, combined with 
the trend between SFR and M», implies a correlation 
between UV luminosity and stellar mass when examined 
over a large range in luminosity. Using the procedure 
outlined above, we recomputed the stellar mass density 
by integrating the UV LF to M AB (1700l) = -18.00 
and allowing the stellar mass to adjust according to 
the empirical relation between SFR and M*. The SFR 
is determined by combining the absolute magnitude of 
galaxies with the luminosity-dependent reddening model. 
The resulting stellar mass densities (Table [6]) suggest 
that roughly as much stellar mass is contained in UV- 
faint galaxies as is contained in UV-bright ones, imply- 
ing a relatively steep low mass slope of the stellar mass 
function, a conclusion that appears to be a generic re- 
sult o f most cosmological simulations (|Nagamine et al.l 
120041 iFinlator etaTI l2007n . Note that this computa- 
tion includes only those galaxies that are brighter than 
M A b(1700A) = -18.00, since it is down to this limit that 
the correlation between SFR and M* has been verified 

8 The zeropoint of the SFR-stellar mass relation appears to 
evolve with redshift ( Nocskc et al. 2007]). For our analysis, we have 
assumed the correlation found at z ~ 2.3. 



empirically. Assuming the correlation is valid at fainter 
magnitudes results in a 64% larger stellar mass density 
contribution from UV-faint galaxies when integrated to 
Mab(1700A) — —16.00. For the subsequent discussion, 
however, we assume the numbers that result from in- 
tegrating to the brighter limit. More importantly, this 
discussion highlights the considerable leeway in adjust- 
ing the stellar mass density estimates upwards even with 
conservative assumptions of the stellar mass distribution 
for UV-faint galaxies. 

9.3.4. Comparisons with the Integrated Star Formation 

History 

Can our revised estimate of the stellar mass density at 
z ~ 2.3 account for the star formation that has occurred 
until then? Figure fT2l shows a compilation of stellar mass 
density estimates from the literature, along with our de- 
terminations and the stellar mass density inferred by in- 
tegrating the luminosity-dependent reddening-corrected 
star formation history. Recall that our calculation of 
the stellar mass density includes galaxies brighter than 
Mab(1700A) = —18.0 at z — 2.3, corresponding to a lu- 
minosity of 0.083L*, where L* is the unobscured charac- 
teristic luminosity at z = 2.3. When integrating the star 
formation history, we must keep track of how these galax- 
ies evolved with time. We have already shown that L* 
evolves strongly with redshift, such that > 0.083L* =2 3 
galaxies at z = 2.3 would have been fainter on average 
at higher z. To account for this fading with increasing 
redshift, we adopt a lower limit to the integral of the 
UV LF that evolves in the same way as L*. Specifically, 
to find the star formation rate density, we integrate the 
UV LF to a limit of 0.083L* where, in the integral, L* 
varies with redshift. By doing this, we are in effect keep- 
ing track of these galaxies' location on the UV LF as 
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a function of time. Note that there may be individual 
galaxies that contribute to the stellar mass density prior 
to z = 2.3, but then fall below a luminosity of 0.083L* at 
some later epoch. However, galaxies will of course accu- 
mulate most of their stellar mass when they are forming 
stars at higher rather than lower rate. Further, virtu- 
ally all of the evolution of the UV LF at z > 2.3 can 
be accommodated by a brightening of L* with increas- 
ing cosmic time (§ I6.1[) and, therefore, it is unlikely that 
there are large numbers of galaxies fading beyond a fixed 
fraction of L* . We noted in § [6] that there is little evo- 
lution in the UV LF between z ~ 4 and z ~ 2.3, but 
presumably most of the galaxies that are fading are des- 
tined to become massive galaxies by z ~ 2.3, and recall 
that our estimate of the stellar mass density at z ~ 2.3 
includes that from massive galaxies with M* > 10 11 M . 
Given these reasons, it is reasonable to assume that in- 
tegrating to 0.083L*(z) should approximate the stellar 
mass accumulated by galaxies brighter than this limit at 
all previous epochs. 

With this premise, we show the integrated star forma- 
tion history in Figure [T2j Accounting for the stellar mass 
content of both UV-bright and faint galaxies results in a 
stellar mass density at z — 2.3 that agrees well with our 
inference from integrating the star formation history. In 
fact, it is perhaps remarkable that with a simple calcu- 
lation where we account for the fading of galaxies with 
increasing cosmic time, we are able to resolve the inte- 
grated star formation history with the global stellar mass 
density at the very epoch where their supposed dispar- 
ity reaches its greatest amplitude (Figure [TO)) . A careful 
analysis of the stellar mass density contributed by galax- 
ies over the bulk of the LF, combined with an integration 
of the luminosity-dependent reddening-corrected star for- 
mation history to an appropriate limit, may obviate the 
need to invoke some other mechanism, such as an evolv- 
ing IMF, to explain the discrepancy. Of course, with the 
present analysis we cannot rule out that there may be 
some redshift evolution of the IMF, and there a re theo- 
retical argu ments as to why this may be the case ([Larson! 
Il998ll2005f ). Indeed, such an evolution may plausibly ex- 
plain the shift in zeropoint of the tre nd between SFR and 
stellar mass as a function of redshift ()Davell2008l ) . All we 
have shown here is that there is a simpler explanation for 
the discrepancy between the integrated star formation 
history and stellar mass measurements at z ~ 2, namely 
that the former must take into account an evolving dust 
correction and the latter are likely to be incomplete for 
galaxies with low stellar masses. In point of fact, incom- 
pleteness of stellar mass density measurements and an 
evolving dust correction are physically well-motivated by 
observations of high redshift galaxies, as we have shown 
here and elsewhere (jReddv et al.ll2006bl R08), whereas 
IMF evolution has yet to be verified observationally. The 
results of the last few sections highlight the subtleties 
that without proper accounting may lead to the types 
of discrepancies reported in the past. Our findings fa- 
vor a more nuanced view of the purported discrepancy 
between the integral of the star formation history and 
stellar mass density measurements. 

Note that we have not measured directly the stellar 
mass function at z ~ 2, but have inferred it by combin- 
ing our knowledge of the UV LF (which gives the number 
density of galaxies) with stellar masses determined from 



broadband SED fitting. Our analysis suggests that an 
appreciable fraction of stellar mass is hosted by sub-L* 
galaxies and that the steepness of the slope of stellar 
mass function may have been underestimated in previ- 
ous studies based on near-IR data. The robustness of 
our conclusions should be verified by significantly deeper 
rest-frame near-IR observations that constrain the low 
mass end of the stellar mass function. A more direct 
sampling of the stellar masses of UV-faint galaxies is re- 
quired. 

9.4. Concluding Remarks 

We conclude this section with a few cautionary re- 
marks. As stated previously, we have adopted the 
z > 4 measurements of the UV LF that are based on 
a maximum-likelihood analysis that is most analogous 
to the method we have used, and that are based on data 
that extend to comparable depths as achieved here, al- 
beit over an area an order of magnitude smaller, in order 
to make the most consistent comparison between star 
formation rate density estimates. Obvious effects that 
can contribute to both random and systematic error in 
the star formation history include cosmic variance and 
the limit to which the UV LF is measured. Large-scale 
multi-field surveys at z > 4 analogous to the present sur- 
vey at z ~ 2 — 3 will provide better constraints on the 
random errors associated with cosmic variance. Further, 
it may be of interest to determine if the similarity in the 
UV LF from the HDF studies versus the universal one 
measured at z ~ 2 — 3 (§ I5.4[) extends to higher redshifts. 

Another point of consideration is the expected turnover 
in the faint-end slope at very faint luminosities. This 
turnover is likely dictated by the threshold of cold gas 
surfa ce density in halos requir ed to trigger star forma- 
tion ()Schmidtlll9"59t iKennicuttl 119980 . In the context of 
the present analysis, the magnitude of the systematic 
effect that results from integrating the LF to zero lumi- 
nosity will depend on the steepness of the faint-end slope. 
For a fixed 4>* , the difference between integrating the LF 
to 0.04L* versus zero luminosity is a factor of 1.85 for 
a = -1.73, 1.44 for a = -1.6, and 1.19 for a = -1.4. 
However, at present there are no empirical constraints 
on the turnover of the faint-end slope of the LF. Even 
locally, where surveys can probe to luminosities signifi- 
cantly fainter than L* , there is no evidence for a fall-off in 
number density of dwarf galaxies. The local it-band LF 
from the SDSS, for instance, appears to abide by a log- 
linear relationship betw een number density and magni- 
tude down to ps 0.02L* (jBaldrv et al.ll2005f ) and, in fact, 
surveys of the local group of star-forming dwarf galaxie s 
suggest an increasing slope to 0.0001L* ([Mated fl998f ). 
Of course, there is no reason why the local rest-frame 
optical LFs should have the same slope as the rest-frame 
UV LF, particularly if the LF represents a sequence in 
mass-to-light ratio. Recall that we have adopted a zero 
lumino sity limit for consiste ncy with the SFRD c ompila- 
tion of iWilkins et~aTl (|2008h and |Hor±kins] (|200l . Given 
the steep faint-end slopes found at z > 2 and the ac- 
commodation of a significant fraction of the luminosity 
density by faint galaxies, we should bear in mind the 
possible systematic effects of integrating to zero lumi- 
nosity, both in terms of the unobscured UV luminosity 
density and the average dust corrections in the case of 
luminosity-dependent reddening. A further systematic in 
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the star formation history may be caused by a redshift- 
dependency in the turnover of the faint-end slope of the 
UV LF. We note, however, that these systematics will 
not affect our comparison of the integrated star forma- 
tion history and stellar mass density measurements at 
z ~ 2 given that we restricted our analysis to luminosi- 
ties (> 0.083L*) where we do have empirical constraints 
on the LF. 

It is also worthwhile to mention that the Kennicutt 
( 1998) relation for converting UV luminosity to star for- 
mation rate is valid only for a stellar population age that 
is > 100 Myr, since it is after this time that the mix of 
O and B stars stabilizes assuming a constant star forma- 
tion history. For galaxies much younger than this, the 
UV luminosity will underpredict the SFR based on this 
relation. Hence, discerning trends in the stellar popula- 
tion age as a function of unobscured UV luminosity is a 
necessary step in computing accurately the star forma- 
tion rate density, particularly since the UV luminosity 
density appears to be dominated by UV-faint galaxies, if 
such UV-faint galaxies are systematically younger than 
their brighter counterparts. 

Finally, we note that because the stellar mass density 
is an integrated quantity, we cannot add arbitrarily large 
amounts of stellar mass at high redshift without violating 
the local constraints, assuming the latter are complete 
in stellar mass. Given that the local measu rements ca n 
be systematically uncertain by up to 30% (|Belll l2003f ). 
and the additional contribution from faint galaxies to 
Mab(17001) = -18.0 at z > 2 is only » 9% of the local 
value, then our finding of significant stellar mass in UV- 
faint galaxies at z > 2 does not pose a problem in terms 
of the budget of stellar mass in the local universe. In 
practice, our preliminary estimates of the number den- 
sity of galaxies with low stellar masses (M* < 10 10 M Q ) 
may be used to constrain cosmological models that cur- 
rently predict low-mass slopes of the stellar mass function 
at z > 3 that are c omparably steep as th e slope of the 
halo mass function (|Nagamine et al.ll2008l ). Ultimately, 
this issue may be resolved through detailed clustering 
analysis of sub-L* galaxies and inferences as to their lo- 
cal descendants. Alternatively, if LAEs represent a short 
phase in the lifetimes of UV-faint galaxies but are other- 
wise unremarkable sub-L* galaxies, then the clustering of 
LAEs may pr ovide clues to the descendants of UV-f aint 
galaxies fe.g.. lKovac et al.ll2007t iGawiser et al.ll2007D . 

10. DISCUSSION: EVOLUTION OF THE FAINT-END SLOPE 

So far, the discussion has focused on what the UV 
LFs can tell us about the star formation rate den- 
sity and buildup of stellar mass. The modulation 
of the LF with respect to the underlying halo mass 
distribution also yields important information regard- 
ing the processes that regulate star formation, such as 
supernovae-driven or radiative winds, and energy in- 
jection from AGN, mechanisms generically referred to 
as "feedback." For example, the sharp cutoff at the 
bright-end of the UV LF may be partly attributable 
to AGN fe edback suppressing star formation in high 
mass halos (ICroton et al.l l 2006l ; IScannapieco et al1l2005t 
iGranato et al.l I2004T ) . even after taking into account 
the saturation of UV light with respect to the to- 
tal star formation rates o f galaxies with larg e SFRs 
(|Adelberger fc Steidell [2005 iReddv et "all l2006bh . Sim- 



ilarly, the shallowness of the faint-end slope of the LF 
relative to that of the halo mass function suggests some 
regulating mechanism associated with star formation it- 
self, such as t hrough reionization (|Kravtsov et al.ir2004t 
iGnedidjgfflp. supernovae winds, or radiativel y-driven 
winds (|Martiijll999t ISpringel fc Hernquistll2003h . 

One important conclusion from our analysis is that the 
faint-end slope of the UV LF is relatively constant and 
steep between z ~ 2 and the highest redshifts where a 
can be measured well, around z ~ 6 (Figure [7]) . At the 
same time, L* evolves strongly between these redshifts 
(§ 16. ip . Because the average galaxy is brightening, the 
invariance of a over these redshifts is not likely reflective 
of an equilibrium between fading and brightening galax- 
ies. Rather, whatever sub-L* galaxies at z ~ 6 brighten 
to become L* galaxies by z ~ 2 are made up in number 
by halos in which gas has newly condensed to form stars 
by z ~ 2. 

The steep a ~ -1.7 at z > 2 stands in contrast 
with the shallower values of a ~ — 1 .1 m easured locally 
(jWvder et al.ll2005t iBudavari et al J 120051) . The redshift 
evolution of a{z) is summarized in Figure[7]and suggests 
that most of the change in a occurs mainly below z ~ 2. 
What may be the cause of this change? It has been sug- 
gested recently that the evolution in a(z), such that a is 
shallower at lower redshifts, may reflect the delayed on - 
set of feedback from Type la SN (jKhochfar et al.ll2007| ). 
However, even if such feedback is energetically impor- 
tant, it is unclear whether it would have any perceivable 
effect on a(z) given that the faint-end population evolves 
strongly between the redshifts in question. A more likely 
explanation is that the evolution in a(z) is dictated sim- 
ply by the availability of low mass halos with cold gas 
at redshifts z < 2. Perhaps it is not surprising that the 
apparent shift from steep to shallow faint-end slopes oc- 
curs at an epoch (z ~ 2) that is marked by a confluence 
of other important transitions, including the reversal in 
the evolution of the cosmic star formation density. 

11. CONCLUSIONS 

We have used the spectroscopic redshifts and pho- 
tometric data in all of the fields of the Lyman Break 
Galaxy (LBG) survey to make the most robust determi- 
nation of the UV luminosity functions (LFs) at 1.9 < 
z < 2.7 and 2.7 < z < 3.4. Our sample includes over 
2000 spectroscopic redshifts, and 31000 LBGs spread 
across 31 spatially-independent fields over a total area of 
3261 arcmin 2 . The depth of these data allow us to select 
LBGs to 0.07L* and 0.1L* at redshift z ~ 2 and 3, re- 
spectively. The LFs are constrained using a maximum- 
likelihood procedure that includes the effects of photo- 
metric errors, contaminants, and perturbation of galaxy 
colors due to Lya. The principle conclusions of this work 
are as follows: 

1. We have quantified the effects of a luminosity depen- 
dent reddening and Lya equivalent width (W-L ya ) distri- 
bution on the incompleteness corrections to our sample. 
Allowing for a larger fraction of galaxies with large Wh ya 
among UV-faint galaxies results in a 3 — 4% increase in 
the faint-end number densities relative to those obtained 
by assuming a luminosity-invariant W^ ya distribution as 
constrained from our spectroscopic sample. Similarly, 
adopting a luminosity-dependent reddening distribution 
where the mean reddening of galaxies decreases to fainter 
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UV magnitudes results in an up to 10% increase in the 
inferred number density of UV-faint galaxies. While 
these differences in the number density are not negligi- 
ble, accounting for these luminosity-dependent system- 
atics does little to alter the Schechter parameters, in 
particular the faint-end slope (a), and it suggests that 
the UV-color criteria are robust to such systematics and 
that our derived LF must be reasonably complete for 
UV-faint galaxies. Adopting reasonable assumptions for 
the luminosity dependence of Wh ya and reddening, we 
derive faint-end slopes of a(z = 2) = —1.73 ± 0.07 and 
a{z = i) = -1.73 ±0.13. 

2. A comparison indicates that our determination 
yields sub-L* number densities that are significantly 
larger, and faint-end slopes that are somewhat steeper, 
than those published previously. We believe our results 
are robust given (a) the large number of spectroscopic 
redshifts used to constrain the bright-end of the UV LF, 
(b) photometry over a large area spread across many 
spatially-independent fields to mitigate cosmic variance, 
and (c) a careful analysis of the systematics (Lya line 
perturbations and luminosity-dependent reddening) that 
are important for computing the faint-end slope. Our 
analysis suggests that LFs based on HDF-N data alone 
are not biased significantly from the universal value of 
the LF determined here, at least at redshifts z ~ 2 — 3. 

3. There is very little evolution in the UV LF in the 
redshift range 1.9 < z < 3.4. However, examined over 
a larger baseline in redshift and using the published re- 
sults at z ~ 6, we find a brightening of the characteristic 
unobscured UV magnitude of ~ 1.2 mag between z ~ 6 
and z ~ 2. The faint-end slope remains relatively con- 
stant and steep between these redshifts, with a value of 
a ~ -1.7 to -1.8. 

4. To examine the frequency of atypical galaxies on the 
faint-end of the UV LF, we compared the number den- 
sity of sub-L* galaxies to those of dusty ultraluminous 
infrared galaxies and galaxies with large stellar masses 
> 10 11 M . With conservative assumptions regarding 
their UV-magnitude distributions, we find that galax- 
ies with large stellar masses and bolometrically-luminous 
galaxies comprise < 2% of the total space density of 
galaxies fainter than 1Z — 25.5. This small fraction un- 
derscores not only the rarity of these objects, but also the 
large number of UV-faint galaxies implied by the steep 
faint-end slope. 

5. Integrating the UV LFs at z ~ 2 — 3 to zero implies 
that 93% of the unobscured luminosity density resides 
in galaxies fainter than L*. Adopting our prescription 
for the luminosity-dependence of reddening, we construct 
bolometric luminosity functions to estimate that > 70% 
of the bolometric luminosity density arises from galax- 
ies fainter than the characteristic bolometric luminosity 
at these redshifts. The luminosity-dependent reddening 
model combined with a steep a imply that the average 
dust corrections needed to recover the bolometric lumi- 
nosity density from the unobscured UV luminosity den- 
sity will depend sensitively on the limit of integration 
used to compute the luminosity density. Of course, these 
corrections will depend also on whether they include only 
the reddening corrections for galaxies routinely selected 
by their rest-UV colors or if they also include corrections 
for galaxies that escape UV selection altogether. 

6. Assuming a constant reddening correction of 4.5 



to the UV-determined star formation history results in 
a factor of two overestimate of star formation rates and 
stellar mass densities accumulated at z ~ 2 — 3 rela- 
tive to the values obtained by assuming a luminosity- 
dependent reddening correction to the star formation his- 
tory. Integrating the latter indicates that at least 25% 
of the present-day stellar mass density was formed in 
sub-ultraluminous galaxies between redshifts z = 3.4 and 
z = 1.9. 

7. The luminosity-dependent reddening-corrected star 
formation history points to a factor of 8— 9 increase in the 
star formation rate density (integrated to zero luminos- 
ity) between z ~ 6 and z ~ 2, significantly steeper than 
the factor of 4 that we would have inferred in the case 
of constant dust reddening. The evolution in the bolo- 
metric star formation rate density is driven equally by 
an evolution in the unobscured characteristic luminosity 
and an evolving (luminosity-dependent) dust correction. 

8. We have examined the offset between the integral 
of the star formation history and previously published 
determinations of the stellar mass densities at z ~ 2, 
the epoch where this discrepancy appears to peak in am- 
plitude and where our data are most sensitive. Given 
the steep faint-end slopes observed at z ~ 2, we have 
explored whether UV-faint galaxies could plausibly ac- 
count for the observed differences. By summing the stel- 
lar mass from all galaxies brighter than 0.083L* =2 i we 
find a stellar mass density that is in remarkable agree- 
ment with the luminosity-dependent reddening-corrected 
star formation history when the latter is integrated to the 
same 0.083L* limit that accounts for the fading of galax- 
ies with increasing cosmic time. This exercise highlights 
the importance of UV-faint galaxies in the total budget 
of stellar mass, and suggests that computing the inte- 
gral of the star formation history in a way that reflects 
how galaxies evolve may obviate the need to invoke other 
mechanisms (e.g., an evolution of the IMF) to reconcile 
the integrated star formation history and the global stel- 
lar mass density at z ~ 2. 

9. Finally, while the faint-end slope at any given red- 
shift is likely to be regulated by feedback, discerning the 
signatures of delayed feedback (e.g., from Type la SN) 
in the redshift evolution of a is not trivial, particularly 
given the strong evolution of the UV LF at z > 2. Our 
results suggest that a is roughly constant at z > 2, con- 
trasting with the shallower values found locally. This 
evolution may be dictated simply by the availability of 
low mass halos capable of supporting star formation at 
z < 2. 
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APPENDIX 
TESTING FOR SYSTEMATIC EFFECTS 

As discussed in § [3J the transitional probability function £ (Eq. [J) is sensitive to changes in the intrinsic properties 
of galaxies. Our goal is to determine if the resulting modulation of £ is significant enough to induce a noticeable 
difference in the maximum-likelihood LF. Since we are interested primarily in computing the UV luminosity distribution 
of galaxies at these redshifts, we will concentrate on how the reddening, N(E[B — V]), and Lya equivalent width, 
iV(WLya), distributions might change with luminosity. More generally, the distributions may also be a function of 
redshift, but for this analysis we ignore such redshift evolution since there is little published evidence for it over the 
redshifts considered here (see below). We conclude by discussing the fraction of stellar objects and galaxies outside 
the redshift ranges of interest. 

The most direct approach for testing systematic changes in N(E[B — V]) and N(WL ya ) is to examine the LBG 
color distribution. However, it is difficult from an analysis of the candidates' colors alone to separate the selection 
effects imposed by the color criteria from those induced by other systematics, such as real changes in the reddening 
and/or W^ ya distributions. A different approach takes advantage of the fact that the transitional probability function 
£ encapsulates all of the information regarding the selection biases imposed by the color criteria. Therefore, rather 
than examine directly the color distribution of BXs and LBGs, we chose to make various assumptions of how the 
reddening and Wi, ya distributions vary as a function of magnitude. Then, based on these assumptions, we recalculated 
£ using Monte Carlo simulations and repeated the maximum-likelihood procedure to find the best-fit LF. 



A. Wiyya Distribution 

There is an increasing body of work that indicates that Lya emitters (LAEs), those objects selected by narrowband 
techniques, exhibit significantly larger rest-frame Wi, ya , but are much fainter in the continuum, on average, than 
traditional color-selected galaxies that are restricted to 1Z < 25.5. Most of these LAEs will lie on the faint-end of the 
UV LF. Adding Lya emission to a star- forming galaxy's spectrum will tend to scatter such a galaxy out of the BX 
selection window (Figure [T3"|). 9 Therefore, for a fixed observed number of faint BX candidates, the incompleteness 
corrections will be larger if they have a distribution skewed towards high W\^ ya , thus increasing the inferred number 
density of faint galaxies. The degree to which the faint-end slope a changes will depend on the W\, ya distribution, 
the stellar population distribution, and the number density of high W-\_ ya systems. For example, there may be little 
difference in the faint-end slope derived assuming a high W\^ ya (LAE) population relative to that derived assuming 
a constant W-^ ya distribution (§ [4]) if LAEs constitute a small and constant percentage of galaxies as a function of 
magnitude on the faint-end of the UV LF. We will now consider in detail how a luminosity-dependent W\^ ya distribution 
affects our analysis. 

A.l. Assumptions 

Here we quantify the effects of a luminosity dependent W\^ ya distribution by making the following assumptions. 
First, for ease of discussion, we assume that an "LAE" is any star-forming galaxy w ith rest-frame Wl vq . > 50 A, 
corresponding roughly to the observational lower limits o f typical LAE surveys (e.g.. lOuchi et al.ll2008h and upper 
< f0% of continuum-selected galaxies to 1Z = 25.5 (e.g., iSteidel et al.ll2003t IShaplev et al.ll2003l . R08). In a strict 
sense, an LAE is any galaxy with Wh ya > 0, but here we limit ourselves to those that are easily identified using 
narrowband techniques, to distinguish them from emission line galaxies that are routinely identified from spectroscopy 
of continuum-selected galaxies. The adoption of the Wi, ya — 50 A cutoff is for reference purposes only, and does not 
affect the subsequent analysis since a separate assumption is made regarding the median value of Wi^ ya for LAEs. In 
particular, the simulations are performed using different values of WL ya rangin g from 50 A t o 250 A (rest-frame), 
the latter being a ca nonical upper limit for standard assumptions of the IMF (Salp eterllI955l ) and solar metallicity 
(ICharlot fc Fal]|[l99l . 

Second, we adopt an average stellar population consistent with the most recent analyses of LAEs at high redshift. 
The range of ages found for LAEs is ~ f Myr at the low end to ~ 1 Gyr at the high end, with typical ages of 
~ 100 — 200 Myr, and low metallicity and reddening (e.g.. lGronwall et "al1l2007t lOuchi et a l. 2008). We conservatively 



9 For galaxies at z > 2.48, where Lya lies in the G-band, the presence of emission will scatter them out of the BX window and into the 
LBG window (see also Figure 4 of R08). 
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Fig. 13. — Color tracks for a typical LBG with constant star formation for 100 Myr and E(B — V) = 0.15 (solid line) and our model for 
Lya emitters (LAEs) with constant star formatio n for 50 Myr a nd E(B — V) = (no reddening; dashed line). The attenuation of colors 
due to the IGM has been accounted for following Madau (1995). The labels along each track indicate particular redshifts as follows: (1) 
z = 1.00, (2) z = 1.68, (3) z = 2.17, (4) z = 2.48, (5) z = 2.65, and (6) z = 2.91. The Lyo line falls in the U n band at 1.68 < z < 2.17 
(between points 2 and 3) and in the G-band at z > 2.48 (above point 4), as indicated by the thicker lines. The effect of adding Ly« emission 
to the spectrum is showp hv-the arrows, lending to scatter galaxies. o,ut_o£ 7 thc BX selection, window J trapezoid) . , r 

assume an average LAE age of 50 Myi " metafheny ol 1 /20 Z© , and zero dust reddening Tor the purposes of our 
simulation. For illustrative purposes, Figure Q1J] shows the colors as a function of redshift for such a stellar population 
compared to a model that represents the typical LBG with age > fOO Myr and E(B — V) ~ O.f 5 (assuming the BC06 
model). As discussed above, irrespective of whether Lya lies in the U n or G band, the effect of adding emission is to 
scatter galaxies at 1.9 < z < 2.7 out of the BX selection window. 

Third, we must det ermine the frequen cy of LAEs among the general star-forming population as a function of 
continuum magnitude. lOuchi et al.l (|2008f ) determine the UV LF of LAEs with W Lya > 50 A at z « 3.1 and find it 



to be fit adequately with a Schechter form of the LF with 0* = (5.6+gJ) x 10~ 4 Mpc~ 3 and M* = = -19.8 ± 0.4 

with a fixed faint-end slope of «lae = —1-6. Integrating the continuum- UV LF of LAEs and comparing with the UV 
LF determined above for all star-forming galaxies implies an LAE fraction (e) that is a strong function of magnitude, 
ranging from less than 0.02% in the brightest magnitude bin (-22.83 < M(1700l) < -22.33) to w 9.3% in the faintest 
bin (—18.33 < Af(1700A) < —17.83). Not e that the fraction of LAEs on the faint-end is based on an extrapolation of 
the UV LF of LAEs assuming a = -1.6 (|Ouchi et al.ll2008D . Further, the fractions will go up (or down) depending 
on whether we decrease (or increase) the limiting Wh ya that segregates LAEs from other galaxies (e.g., using a limit 
of W Lya — 20 A instead of 50 A). In our simulations, we assume (1) no evolution in the UV LF of LAEs between 
z ~ 2 — 3 and (2) a fraction of LAEs that varies with absolute magnitude in accordance with our findings above, with 
a fraction of 10% for Wh ya > 50 A in the faintest bin considered here. 

To recap, the main assumptions going forward are that LAEs are described by a 50 Myr stellar population with no 
reddening and comprise anywhere from < 0.02% to 10%, respectively, of galaxies within the bright and faint magnitude 
bins of our analysis. In the next section, we consider how different values of W-L ya among LAEs affects the faint-end 
of the UV LF of all star- forming galaxies. 

A. 2. Effect of a Changing Whya Distribution 

With the aforementioned premises, the UV LF is computed for varying amounts of emission among the LAEs, with 
equivalent widths from 50 A to 250 A. In our simple model, the luminosity dependent W Lya distribution can be 
expressed as: 

N(W Lya , U) = [1 - e(R)\ N a (W Lya ) + e(K)5(u), (Al) 

where ^(lULya) is the distribution for 1Z < 25.5 galaxies (Figure[T]), e(lZ) is the WL yQ > 50 A fraction as a function 
of magnitude as determined above, and S(lo) is a delta function with center at u> = 50, 100, 150, 200, and 250 A. 10 
The results at z ~ 2 and z ~ 3 (Figure [14]) are presented in terms of the ratio of the number density (77*) for 
different values of W Lya for LAEs to the number density (77) derived using the fiducial W Lya distribution for 1Z < 25.5 
continuum-selected galaxies (e.g., Figured]), 



'L 
v 



(A2) 



10 The intrinsic distribution for UV-bright galaxies (Figure [TJ also includes a small fraction of continuum-selected galaxies with WL ya > 
50 A. We ignore this small overlap between the continuum and LAE W-L ya distributions, with the obvious consequence of slightly increasing 
the total fraction of galaxies with Whya > 50 A. 
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Fig. 14. — Comparison of maximum-likelihood number density of galaxies as a function of magnitude for (1) our fiducial UV LFs assuming 
that all galaxies abide by the W^yo distribution seen for bright (TZ < 25.5) continuum-selected galaxies (Figure [TJ and (2) the UV LFs 
derived assuming a population of Lyce emitters (LAEs) with high WL ya at faint magnitudes. The different lines show the ratio of the LFs 
determined from (2) to that determined from (1), and correspond to differing amounts of Lya emission attributed to the LAE population. 
The typical error in this ratio (<Jf) is shown by the vertical errorbar. In all cases, we find that such a population of high WL ya systems 
does little to alter the faint-end slope of the UV LF. 

with error given by 



2^1/2 



(A3) 



The error in number density cr^« is determined in exactly the same manner as the error in the LF (c^; §H]). Even with 
the most conservative assumption for the Whya attributed to LAEs (an assumed value of Whya = 250 A), we find 
that inclusion of such a population alters the faint-end number densities at a 2 — 4% level depending on the redshift. 
Not surprisingly, at lower redshifts, 1.9 < z < 2.7, the inclusion of LAEs increases the inferred number densities at 
the faint-end by < 3%, since such galaxies would be preferentially scattered out of the BX selection window. The 
opposite is found at higher redshifts (2.7 < z < 3.4), where the number density is systematically lower by up to 4%, 
primarily because there are more LAEs scattered from z < 2.7 into the LBG selection window than are (at any redshift 
2.7 < z < 3.4) scattered out of the window, assuming no evolution of the UV LF of LAEs at the redshifts of concern. 
Hence, for a fixed number of z ~ 3 LBG candidates, the tendency would be to over-estimate the number density had 
we not accounted for the LAE population. 

The expectation of a steepening a for the universal UV LF (i.e., for all star- forming galaxies) at z ~ 2 when including 
LAEs is not borne out with our simulations for several reasons. First, the overall fraction of LAEs, even at the faint- 
end of the UV LF, is small (< 10%). Their effect on the LF is further diminished because they will only affect the 
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TABLE 7 

Best-fit Schechter Parameters for UV LFs of 1.9 < z < 3.4 Galaxies 



Redshift Range Systematic Effect a M^ B (1700A) 0* (xlCT 3 Mpc -3 ) 



1.9 < z 


< 2.7 


Fiducial a 
W hya = 150 A b 

[E(B - u)] grad c 

Combined 1 ' 


-1.67 ±0.06 
-1.68 ± 0.06 
-1.72 ± 0.06 
-1.73 ± 0.07 


-20.65 ±0.08 
-20.68 ±0.08 
-20.71 ±0.09 
-20.70 ±0.11 


2.84 ±0.42 
2.83 ±0.42 
2.64 ± 0.46) 
2.75 ± 0.54 


2.7 < z 


< 3.4 


Fiducial a 
W hya = 150 A b 

[E(B - v)] grad c 

Combined d 


-1.77 ±0.12 
-1.80 ±0.12 
-1.78 ± 0.12 
-1.73 ± 0.13 


-20.98 ±0.13 
-21.06 ±0.13 
-20.97 ±0.13 
-20.97 ±0.14 


1.54 ±0.43 
1.34 ±0.42 

1.55 ± 0.43 
1.71 ± 0.53 



a Fiducial LF assumes no change in the WL ya and E(B — V) distributions as a function of 
UV apparent magnitude. 

b LF derived assuming an LAE population at the faint-end with WL ya = 150 A (see text). 
c LF derived assuming a linearly declining mean E(B — V) for galaxies with 1Z > 25.5 (see 
text). 

d LF derived combining the effects of a changing WL yQ and E(B — V) distribution as a 
function of UV apparent magnitude. 



broadband colors if they are at partic ular redshifts (F igure [T3")) . Second, the current determination of the UV LF of 
LAEs assumed a fixed «lae = —1-6 (|Ouchi et al.1 120081 ) . which is not too different from the fiducial a (§ QJ. Thus, 
while the fraction of LAEs among the general population is a strong function of magnitude between our brightest and 
faintest magnitude bin, it is in fact a relatively constant « 8 — 9% for bins fainter than M(1700A) k, —20. Adopting 
a steeper faint-end slope of the UV LF of LAEs of Q!laelf = — 1.9 still results in a universal faint-end slope (auvLF) 
that is consistent with the fiducial value (Figure [T4f . 

Future studies that constrain more robustly «laelf over the entire redshift range probed here will be useful for 
assessing the overall impact of high W^ya systems on the faint-end slope inferred for all star-forming galaxies. It is also 
not unreasonable to suspect disparate absorption properties between UV bright versus faint galaxies, so spectroscopic 
studies, while diffic ult to carry out at present, are crucial for assessing the variation of W^ya with luminosity (e.g., 
iShaplev et afll2003h . 

Nonetheless, while the inclusion of high W\^ ya systems among the UV continuum-faint population may be a small 
systematic effect (3 — 4%), it is not negligible compared to the error in number density at the faint-end of the UV LF 
(10 — 15%), and so should be included in any proper assessment of the UV LF. The critical point, and one that is 
demonstrated unambiguously with our simulations (Figure [T4"]) . is that such a systematic effect does little to alter the 
faint-end slope of the universal LF. 11 In our final determination of the LF, we have made the conservative assumption 
of a median value of Wh ya — 150 A for the LAE population with Wl Y q > 50 A (although, as noted above, the exact 
value does little to alter the LF). The resulting best-fit Schechter parameters are listed in Table [71 

B. REDDENING DISTRIBUTION 
B.l. Test Cases 

In the prior section, we assumed a young stellar population and no reddening when modeling LAEs. In this section, 
we test for modulation in the LF if all UV-faint galaxies are characterized with a young stellar population and lower 
reddening than UV-bright galaxies. In this case, the LAEs would simply represent a phase of UV-faint galaxies with a 
short duty cycle of < 10%, based o n the number density o f LAE s compared to the g eneral continuum population (e.g., 
iKovac et al.Tl2007t iNagamine et al.ir2008t IVerhamme et alJl2008HGawiser et al.ll2007l ). Note that here we are concerned 
with the changing distribution of reddening of galaxies as a function UV magnitude at a given epoch. A somewhat 
related issue is how the reddening distribution in general shifts to lower values at higher redshift for galaxies of a given 
star formation rate (R08; see also §[8]). 

The first case under consideration is if the average E(B — V) of galaxies with 7Z > 25.5 is zero. This scenario would be 
the most conservative one we can make for two reasons. First, R08 demonstrate using UV continuum-slopes and Spitzer 
MIPS data that the reddening distribution of UV-selected galaxies does not vary significantly to the spectroscopic 
limit of 1Z ~ 25.5. Because this limit is arbitrary with respect to galaxy properties, wc would not expect the reddening 
to change so suddenly for galaxies fainter than this limit and indeed it would be unphysical. A more meaningful model 
is one in which the average reddening asymptotes to zero proceeding to fainter luminosities (R08). Second, there is a 
non-negligi ble fraction of galaxies at z ~ 2 — 3 on the fain t-end of the UV LF that are bolometrically-luminous and 
dusty (e.g.. lChapman et al.ll2005l : Ivan Dokkum et al.ll2004h . Yet, a large fraction o f these galaxies have colors and UV 
opacity that do not differ significantly from those of UV continuum-bright objects (|Chapman et al. 2005: Rcdd v~t al.l 
2005. 12006U R08). By assuming an average reddening that falls to zero for galaxies fainter than 1Z = 25.5, we are 
effectively seeking the maximal change in the LF under such a scenario. We also consider a more physical reddening 

11 While the small change in the W\^ ya distribution brought on by the inclusion of LAEs does little to alter the faint-end slope, significant 
discrepancies in the faint-end of the UV LF arise when not accounting at all for the W\^ ya distribution of galaxies at these redshifts (R08). 
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distribution whose mean (E(B — V)) decreases monotonically from a value of (E(B — V)) ~ 0.13 at 7Z = 25.5 to 
zero at the faintest apparent magnitude bin of our analysis (case 2). With this model, galaxies at the faint-end of the 
UV LF (M(1700A) fss —18.00) will have close to zero reddening, similar to the mean redden ing (as inferred from the 
rest-UV slopes) of dropout galaxies at z ~ 6 with comparable unobscured UV luminosities (|Bouwens et al.ll2006f h If 
we define the N(E(B — V)) distribution for galaxies brighter than 1Z = 25.5 as 

N[E(B -V),K< 25.5] = N a , (Bl) 

then our model for the luminosity dependence of N(E(B — V),7Z) can be written as 

N(E(B -V),K)=N , ft < 25.5 

= f((E(B-V)),a(N )), 
ft > 25.5, (B2) 

where the function f{(E(B - V)),cr(N )) is a Gaussian with mean (E(B - V)) = (Case 1) and (E(B - V)) = 
—0.097?.+ 2.43 (Case 2) and dispersion equivalent to that observed for 7V (i.e., a(N Q )). Note that while N Q is not in 
fact distributed normally, the differences that arise by assuming a Gaussian are negligible. Further, for simplicity we 
have assumed that the dispersion of the E(B — V) distribution is independent of magnitude. R08 argue for an increased 
dispersion at faint magnitudes due to the mixing of galaxies with intrinsically low star formation rates and those that 
are UV-faint because of high extinction. The effect of such an increasing dispersion is to reduce the effective volume 
of the survey and thus the incompleteness corrections will be larger at the faint-end. In general, however, because the 
number density of UV-faint galaxies with little reddening is inferred to be much larger than that of heavily reddened 
UV-faint galaxies (§[7]), the increase in dispersion attributable to the latter is likely to be negligible. In addition, at the 
faintest magnitudes where {E(B — V)) approaches zero, the dispersion will be dominated not by reddening but by the 
intrinsic variation in SEDs of galaxies. Consequentially, the dispersion will likely be smaller than the <j(Nq) measured 
at brighter magnitudes. Note that there is a non-negligible number of galaxies in the spectroscopic sample that have 
measured E(B — V) < 0. Since we use E[B — V) as an indicator of dust, we set a lower limit of E{B — V) = for 
any galaxies that happen to be assigned a negative value. Finally, R08 demonstrate that the mean extinction among 
galaxies above a particular unobscured luminosity is roughly constant with redshift between 1.9 < z < 3.4. Motivated 
by this, we adopt non-evolution of reddening in the simulations; i.e., N(E(B — V), ft, z) » N(E(B—V),1Z). For brevity 
in the subsequent discussion, the abbreviation ZR refers to the case of a discontinuous reddening distribution such that 
galaxies with ft > 25.5 have E(B — V) = 0. Similarly, LDR refers to our analytical model for the luminosity-dependent 
reddening distribution that has reddening decreasing monotonically with UV magnitude. 

B.2. Results 

The results are summarized in Figure [15] for both case 1 (ZR), where all galaxies with 1Z > 25.5 have E(B — V) = 0, 
and case 2 (LDR), where the mean reddening decreases with magnitude. There are several conclusions of import. 
Focusing on the lower redshift galaxies, in the ZR case we find a significant increase of up to 70% in the inferred 
number density of galaxies with M(1700A) > —20 relative to that inferred from the fiducial model. This can be 
understood by examining Figure [T3"l Galaxies with no reddening (e.g., as in the case of the LAE model) have bluer 
colors that approach the boundary of the BX color selection window. Thus, photometric errors preferentially scatter 
galaxies out of the selection window compared to typical BXs. The net effect is that for a fixed number of observed BX 
candidates, the incompleteness corrections will be larger, leading to larger number densities. We find a statistically 
insignificant difference between the fiducial and ZR faint-end slopes, attributable to the fact that the reddening 
distribution is fixed to have (E(B — V)) — 0. A gradually declining distribution (case 2; LDR), results in number 
densities that are up to « 10% larger and a slightly steeper faint-end slope, although the parameters of the Schechter 
function are still consistent with those of the fiducial case within their respective marginalized errors. 

The luminosity-dependent variation of the E(B — V) distribution has less of an effect on the faint-end number 
densities in the higher redshift range 2.7 < z < 3.4, primarily because the LBG criteria include colors that are much 
bluer than those expected for even a young and unreddened stellar population. However, because we account also for 
the intrinsic scatter in N(E(B — V)), the differences in faint-end number densities are still somewhat larger than we 
would have obtained had we modeled the N(E(B — V)) distribution as a 8 function at a given ft-band magnitude 
(i.e., assuming all galaxies at a given 1Z have the same reddening). In any case, the LBG criteria are somewhat more 
robust to extreme assumptions regarding the E(B — V) distribution of UV-faint galaxies compared to the BX criteria. 
Simply altering the BX criteria to include candidates with bluer G — 1Z colors could alleviate some of this difference, 
but we note that our deep photometry indicates that galaxies with such blue colors (G — 1Z < —0.2) are rare. For the 
purposes of our present analysis, the exact placement of the color criteria is not important as long as the effect of the 
criteria is modeled appropriately and incompleteness is accounted for using a likelihood analysis. 

Deep spectroscopy combined with multi-wavelength indicators of the extinction of sub-L* galaxies is required to 
more robustly constrain the E(B — V) distribution as a function of unobscured UV luminosity. For the time being, 
however, we have shown that adopting simple scenarios for how the distribution changes with UV luminosity results 
in a faint-end slope that can be potentially steeper at z ~ 2 than we would have obtained by assuming a constant 
E(B — V) distribution. For our final determination of the UV LF, we have adopted our prescription for the luminosity 
dependence of reddening, namely one in which the reddening declines linearly with apparent magnitude, as discussed 
above. The resulting best-fit Schechter parameters are listed in Table [7l 



28 



Reddy & Steidel 



i i | i i i i | i i i i | i i i i | i r 




1Q -4 I I I I I I I I I I I I I I I I I I I I I I 1Q -4 I I I I I I I I I I I I I I I I I I I I 

-21 -20 -19 -18 -21 -20 -19 -18 

M JB (1700A) M aB (1700A) 



Fig. 15. — Change in faint-end slope of the UV luminosity function assuming that galaxies with TZ > 25.5 have (1) {E(B — V)) = 
( a E(B— V)=o) an d (2) (E(B — V)) that falls off linearly with magnitude (<*Sg_ v pi compared to the fiducial value that assumes galaxies 
have the same (E(B — V)) ~ 0.13 irrespective of apparent magnitude («o). For clarity, data points are shown only for case (2) in both 
panels and the Schechter fit for case (2) is excluded from the right panel. 



C. Objects Outside the Redshift Ranges of Interest 

Intrinsic variations in the SEDs of star-forming galaxies and photometric errors lead to significant wings of the 
redshift selection functions, N(z, L), for color-selected samples. Here, the selection function is parameterized in terms 
of both redshift and luminosity. In this analysis, we have computed the LFs specifically over the redshift ranges 
1.9 < z < 2.7 and 2.7 < z < 3.4, although the exact redshift interval used is unimportant since the LF does not evolve 
over these redshifts (see R08 and §[6]). More importantly, we estimate the fraction of photometric candidates that fall 
in the redshift ranges over which the LFs are computed, information that is provided directly from the spectroscopic 
sample. 

In the faintest apparent magnitude bin of the spectroscopic sample, 25.0 < TZ < 25.5, the observed fractions of 
galaxies (excluding AGN/QSOs) that fall in the ranges 1.9 < z < 2.7 and 2.7 < z < 3.4 are 77% and 72%, respectively 
(virtually all of the objects that are outside these redshift ranges still have z > 1 since the contamination due to z < 1 
objects is very small at these faint magnitudes — see R08). In the previous analysis, the fractions of 77% and 72% are 
assumed to remain constant for galaxies fainter than TZ — 25.5. There are at least a couple of reasons why this is likely 
to be a reasonable assumption. First, the contamination from z < 1 objects to the photometric sample is a strong 
function of magnitude. If this trend continues to fainter magnitudes, then we would expect the z < 1 contamination 
rate to be less than 1% for objects fainter than TZ = 25.5. Adopting zero contamination from z < 1 objects at TZ > 25.5 
results in number densities at the faint-end that are f=a 1% larger. This is an insignificant change given the magnitude 
of the other systematics discussed above. In theory, the larger photometric errors for UV-faint objects may result 
in an increase of the contamination rate at the faint-end . A sm all sample of 15 spectroscopically-confirmed galaxies 
with 25.0 < TZ < 26.0 in the Q1422 field (jSteidel et all [2003; R08) implies a zero contamination fraction. Larger 
spectroscopic samples of UV-faint objects will be required to obtain a statistically-robust estimate of contamination 
at the faint-end. However, all of the previous s tudies that have attempted to constrain the faint-end of the UV LF 
(jSteidel et al.l Il999fc ISawicki fc Thompson! I2006D also assume negligible contamination at the faint-end. Hence, the 
difference between our determination of a steep-faint end slope and the shallower values found elsewhere (see § 15. ip 
cannot be attributable to differences in the assumed contamination rate as a function of UV luminosity. 

Second, the redshift distribution for UV-selected galaxies will be modulated primarily by systematics that affect 
the overall colors of galaxies at these redshifts, namely Lya perturbations and the E(B — V) distribution, such that 
N(z,L) = f(N[Wi iya ],N[E(B — V)]). From this discussion, we conclude that the redshift distribution likely remains 
similar between UV-bright and faint galaxies (i.e., N(z,L) « N(z)). It is easy to see, however, the potential danger 
of assuming that N(z) is insensitive to luminosity: even gradual changes in the stellar population and reddening of 
galaxies as a function of magnitude may result in an artificial suppression of the faint-end of the UV LF with respect 
to the bright-end. This motivates the need for selection criteria that are efficient at targeting galaxies with a wide 
range of intrinsic properties at the desired redshifts (see next section). 

For the selection criteria adopted here, the presence of high W^ ya systems and/or a bluer population of UV-faint 
galaxies does little to alter the parameterization of the maximum-likelihood LF under reasonable assumptions for LAE 
number densities and the reddening distribution. In other words, the modulation of N(z,L), and more generally £, 
the transitional probability function, due to these systematic effects do not affect appreciably our LF determination. 
Substantial changes in the redshift distribution can arise from a rapid evolution of the LAE number density and 
E(B — V) distributions over the redshifts of concern. However, lacking evidence that such evolution is occurring — 
and, as discussed above, there is little evolution in the E(B — V) distribution over these redshifts; R08) — it is likely 
that the redshift distributions of BXs and LBGs with TZ > 25.5 is similar to those of TZ < 25.5 galaxies. 
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D. Implications for Color Selection at High Redshift 

It is instructive to take a broader view and examine the significance of the tests described here in the context of 
the color selection criteria. The primary result of this section is that the systematics brought about by reasonable 
assumptions for the unobscured UV luminosity dependence of the W\^ a and reddening distributions do little to alter 
our inference of the UV LF for the selection criteria used here. We emphasize the latter part since obviously some sets 
of criterio n will be more susceptible to modulations of N(W i, V a) and N(E[B — V]) than others. As has been discussed 
elsewhere (Steidel ct al. 2003, 2004; Adclbcrgcr et al. 2004), the goal of observing efficiency dictates that a balance be 
struck between the inclusion of as many galaxies at the redshifts of interest as possible and the exclusion of as many 
foreground or background contaminants. In the context of the present analysis, luminosity-dependent properties of 
galaxies should also be taken into account when designing color criteria. UV-dropout criteria are in general the most 
efficient method for selecting high-redshift galaxies. Because they target preferentially bluer galaxies, the luminosity- 
dependent systematics expected for UV-faint galaxies works in favor of their selection via rest-frame UV emission. In 
contrast, near-IR selections that target redder galaxies (either because they are dusty, have large stellar masses, or 
both), may potentially miss an appreciable fraction of galaxies that populate the faint-end of the UV LF (§[7]). Hence, 
the aggregate of these selection methods provides a complementary view and are necessary for obtaining an unbiased 
census of star formation. 

Obviously, incompleteness corrections for criteria that target blue galaxies are not completely immune to extreme 
luminosity-dependent gradients in N(E[B — V],1Z), for example (Figure Il5j) . The power of our simulations and 
maximum- likelihood method is that they can be used to quantify and correct for even severe biases (e.g., in the faint- 
end slope) imposed by high redshift galaxy selection in a way that is not possible with the traditional methods of 
computing luminosity functions (see discussion in R08). The keystone of our entire method is the spectroscopy: while 
not extending fainter than the typical ground-based magnitude limit, spectroscopy for UV- bright galaxies does provide 
a critical foundation, or "zero-point," for assessing how luminosity dependent systematics may affect our inferences of 
the faint-end. Using these techniques, we showed in R08 that, after correcting for low redshift objects and AGN/QSOs 
based on extensive spectroscopy, the UV LF inferred by magnitude limited surveys is similar to that derived from 
color-selected surveys. Hence, we argued that we must be complete for UV-bright galaxies. The slight modifications 
of the Schechter parameterization of the LF after taking into account N(W^ ya ,TZ) and N(E[B — V],1Z) implies that 
our determination must also be reasonably complete for UV-faint galaxies. 



